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1.0 Introduction 



The problem to be considered is that of approximating functions 
of several variables. We will assume that values of the functions are 
known on a cartesian grid of points, say at / , . . . ,X? \ 

V ^2 \)* 

ij = l,...,Nj, j = l,...,n. (Note: We will often use superscripts, 

which should not be interpreted as exponents) . We will denote the function 

values by F . . 

1 n 

The application in which we are particularly interested is the 
representation of aerodynamic and infrared signature data for use in 
real time hybrid computer simulations. Thus the approximations must be 
capable of being evaluated rapidly. There are usually large amounts of 
data available and it is desirable and sometimes necessary to reduce the 
amount of data which must be stored. 

The data is obtained by measurement, hence may be subject to error, 
in some instances as much as 10% or more of the range of the function. 

This fact will have its impact on some of the measures of the goodness of 
fit which we seriously consider. 

There has been theoretical work done on this problem, usually in a 
more general setting. Notably there is the work of Diliberto and Strauss 
[1] and Golomb [4] , as well as some more general and some less general, such 
as [2], [5], and [7]. There are also a number of practical articles in the 
area which have had an influence, such as [3] and [6]. 

We will briefly discuss several common methods using the full data 
set. We will be more interested in approximation schemes which reduce 
the amount of data and will discuss some of these in more detail. Some 
Fortran IV programs for processing the data into the appropriate form 
are given in the appendix. 
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2-0 Interpolation on the Full Data Set 
2-1 Linear Interpolation 



The simplest and most obvious way to approximate the function 
using the full data set is to use iterated linear interpolation 
between points. This results in the approximating function being 
multilinear (i.e., linear in each variable). This is sometimes called 
characteristic strip interpolation. The approximation is continuous, 
but does not have continuous derivatives in general. 

A simpler version of the above uses a hyperplane to fit the function 
and its first divided difference in each variable. This is equivalent 
to a linear Taylor series representation where the differences replace 
the derivatives. This approximation is easier to compute, but it is not 
continuous. 

The two above approximations have been available in hardware from the 
MacNeal - Schwendler Corporation [8]. A digital processor is used to com- 
pute a function of up to four variables. The input variables are analog 
and are converted to digital. The output function is available in both 
digital and analog form. 

It should be noted that the hyperplane fit can be arranged so as to 
give a continuous approximation. Instead of using the same hyperplane 
over the entire n-cube, one uses many. After the n-cube in which the 
point lies has been determined, one then must determine a subset of n 4- 1 
vertices which define the n-simplex in which the point lies. (It will be 
unique, unless the point lies on a boundary of one of the simplexes) . The 
n + 1 values of the functions at the vertices of the simplex then define 
the hyperplane. This procedure is also available in hardware from the 

MacNeal- Schwendler Corporation, although the author has no other information 
about it. 
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2. 2 Smooth Interpolation 



For some applications it may be necessary to provide approxi- 
mations which are smooth. One could then consider some type of 
global spline approximation. Because of the necessity of rapid 
evaluation it would be necessary to preprocess the data. This would 
result in an increase in the amount of stored data and in many cases 
this would not be feasible. One could also use a local approximation, 
such as Hermite interpolation, or deficient quintic splines . For 
rapid evaluation these would also need to be preprocessed, with an 
increase in the amount of stored data. 

Because of the increase in the amount of data it is probably not 
feasible to use the full data set when smooth approximations are desired. 
If the approximation is to be smooth in one or two variables, the in- 
crease in the amount of data may be manageable. The procedures in these 
cases are well documented and will not be discussed here. 

3.0 Approximations 

The type of approximations we will consider here are generally of 
the form of sums of functions of fewer variables. This can be a very 
effective means of fitting certain functions, and certainly is a de- 
sirable method from the standpoints of both evaluation time and storage 
requirements. However, for it to be an accurate method the function to 
be approximated must have some special properties. In particular, it 
is necessary for the function to have basically the same character for 
different cross sections. For example, consider the approximation of 
F(x,y) by f (x) + g(y)* For fixed y^ the approximation becomes 
F(x,yQ) - f (x) + g(yQ). If F(x,y) has basically the same shape for 
various constant y, then the approximation will be a good one. 
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If F(x,Yq) and F(x,y^) have basically different shapes, then the 

approximation will not be good. We will discuss several schemes for 

obtaining approximations by sums of functions of fewer variables. 

There are many possible variations on this theme. For example, 

one could approximate a function of four variables by sums of products 

11 2 2 

of functions of two variables, e.g. f (w,x) g (y,z) + f (w,y) g (x,z) 

3 3 

+f (w,z) g (x,y). This replaces a function of four variables by six 
functions of two variables. Unless one of these products can be 
eliminated, from the standpoint of evaluation time one is not likely to 
gain much, although the storage required will likely be much less. 
Nonetheless, the above method appears like it could be fruitful, and 
should be explored. The form of the approximation has the advantage that 
for functions of three variables it reduces to a sum of products of func- 
tions of two variables and ones of one variable. 

When the approximations are evaluated, it is now much easier to use 
smooth approximations, since the number of variables has been decreased. 
3.1 Minimax Approximation by Sums of Functions 

12 n, 

We consider the approximation of a function F(x ,x , . . . ,x ) by 

^f^(x^) and seek to minimize , sup |f(x^, x^,...,x^) - ^f^(x^)| over 

i 

12 n 

all possible functions f , f ,...,f . In our case we have a discrete 
function, so we wish to find vectors ^f^,...,f^ ^^ . . . ,|f^, . . . ,f^^^ 



which minimize 



^j.) ' 



max 



|f. , . I over all possible 

12 n j j 



An algorithm for the solution of this problem was given by Diliberto 
and Strauss [1] and it was analyzed, and convergence proven by Golomb [4]. 
A sequence of functions 



(vectors) , ( » • * • 
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£ = 0,1, , is constructed ao that |f|, . . . ,f^ ) ) 

is a solution to the problem. In general the solution is not unique. 

However, the above solution has the desirable property the maximtim error 

for any fixed index is minimized. The algorithm can also be applied if 

one or more of the X*? take on their values in multidimensional space. 

Hence, the algorithm could be used to calculate the best approximation 

to a function of six variables, F, by a sum of the form 

i- i^ 1. le- ^ 

, 0 ZiSA 123t-56 

C + ff’ ; + £. . : , for example. 

If i2'^3 „ 

The algorithm is as follows. Initialize f^* = 0, i. = 1,...,N., 

0 j ^ ^ 

i = and F, = F, for all . Then 

’ i, 1 i ’ 1 12 n 

Inin 

for k = 0, l,...,set j = (k mod n) + 1 and Z = [k/n] + 1, where denotes 



the integer part. 


For i . = 
J 


1 , 2 , . . . ,N , let 
J 




Af = h 


l^max 


F^ 

F.i 

12 n 


min 




l£ i £N 




Ki <N 

— n — 1 



1 n 



1^ p_<n 

P^j 



1 <_ pj<n 
P^j 



k+1 k 

h ^n ^1 ^n - p - p* -P- . P J 

fj’"^ = + Af. 

j j 

Although not analyzed by Golomb, a similar algorithm for minimax 

12 13 

approximation by sums of functions of two variables, f.* ^ f ^ ^ ^ 

, , ^1’ 2 1’ 3 



+ + 

i 1 
^1’ n 



+ f^ l>n been programmed and successfully 

n-1 n 



applied to many data sets. This algorithm is implemented in program 
MINMAX, whose use is described in the appendix. 

3. 2 Least Squares Approximation by Sums of Function 

We now turn to least squares approximations. As in the previous 
instance we will consider fitting with sums of arbitrary functions of 
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fewer variables. The easiest form of least squares approximation is 



an analysis of variance technique outlined by Pike and Silverberg [6]. 

The algorithm is quite simple to apply. One sequentially fits functions 

of zero (a constant), one, two, etc., variables. We proceed as follows. 
^ N, N 

Let K = 



n 



HN. i =1 



1 n 

I •••If 



, the mean value, and let the 



i =1 



n 



n 



residual R. > 
^1 



= F. 
n 1 



- K, 1 < i. < N., 1 < 1 < n. 

1 --J- 



n 



For each j = we define 

N, N 



tJ = 



n 



1 j-1 

I ... I I 



N 



N... 

^ ^ r"" 1 

... R , , . • . , , for 1 , 1 . 

i =1 i. 1=1 1 •4.1=1 i =1 1 ^ 

UN, 1 j-1 j+1 n 



k=l 



• N, 



That is, is the average value of r|^ over all values with j— 

subscript ^ fixed at i.. ^ ^ 

J 



2 1 / 1 2 
Then the new residual is R. =R.,...,. “(V. +V. 

1. 1 1 1 I 1 1, 

1 n 1 n ' 1 z 



n 



+ ... + V . 

1 



A sum of all possible functions of two variables is fit in exactly the 

same fashion, taking the value for a particular pair of subscripts to be 
2 

the average of R. over all values with that pair of subscripts 
1 n 

fixed. 

One desirable feature of this procedure is that the sum of the 
variances of the fitting functions and the residual is equal to the 
variance of the function itself. For example 



n 




Thus one can easily see how much each component is contributing to the 
approximation, Those which do not make a significant contribution may be 
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Unfortunately there seems to be no general rule for observing data 
and then selecting an appropriate approximation scheme. While plots of 
the data may help, there is sometimes not a satisfactory answer. 

Tables 1-6 give the data on which the programs were tested. Tables 
7-12 give the pertinent information about the approximation: Type of 

approximation, maximum error, root mean square error (for least squares 
approximations), and the approximate evaluation time in terms of the 
time required for a one dimensional table look up. In determining the 
last number, it was assumed that an n-^dimensional table look up takes 
2^-1 one dimensional table look ups. In cases where more than one 
table is used, some time will be saved since the argument processing 
needs to be done only once. 

Data Set 1: This data represents the function 

f(w,x,y,z) = exp^(w+l)^ - x/(l+y)j on a uniform grid of 5x9x11x11 points 
on the unit 4-cube. This data was selected to test the programs on a set 
of smooth data which was not a sum of functions of fewer variables, and 
which behaved in a different but well behaved manner in each variable. 
Inspection of Table 7 shows that a quadratic fit in any of the first 
three variables results in a good approximation. None of the sum of 
functions of fewer variables is good. Use of a quadratic fit in the 
third variable results in the data being reduced by over 70%, with a 
possible increase in evaluation time, compared with four dimensional table 
look ups, depending on the argument processing time. 

Data Set 2: This data set represents actual test data obtained from NMC. 

Cross-section plots for the various constant values of the second and third 
variables were available. Perusal of these revealed that the variation 
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in the first variable would probably be fit nicely by a quadratic. 

However, the cross-section did not differ by only a constant, thus it 
appeared a sum of functions of one variable would not give a satisfactory 
approximation. This was indeed the case, as Table 8 shows. The smallest 
rms error was achieved by a quadratic fit in the second variable. The 
quadratic fit in the first variable was almost the same, and resulted in 
the amount of data being reduced by over 85%. Evaluation might be slightly 
longer than 3 dimensional table look up, depending on argument processing 
time. For this data set only we also tried a linear approximation. This 
revealed that the smallest rms error was for the fit in the second variable. 
This would result in the amount of data being reduced by over 75% and a 
slight reduction in evaluation time compared with three dimensional table 
look up. It is interesting to note that in the third variable the maximum 
error is smaller for a linear fit than it is for the quadratic fit. 

If one wished to consider a sum of two variables approximation, it is 
observed that the maximum error obtained by MINMAX was .164, while for 
APPROX it was .219. Near the center of the table the errors are basically 
of similar size, although MINMAX gives the minimum possible error over any 
two fixed indices in this case. 

Data Set 3: This data set is very similar to the previous and the same 

comments apply. See Table 9 for the details. 

Data Set 4; This data set also represents test data obtained from NMC. 
Cross-section plots revealed it to be unlikely that any of the approximations 
will be very good. The best is the quadratic approximation in the second 
variable, which reduces the amount of data by two- thirds, but increases the 
evaluation time slightly. None of the sum of functions approximation are 
worth considering. See Table 10 for the details. 
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Data Set 5: The comments on this data set are basically the same as 



on the previous set. See Table 11 for the details. 

Data Set 6: This data also represents test data obtained from NMC. In- 

spection of the data reveals that a sum of functions approximation will 
probably not be very good. Table 12 verifies this. The best is a 
quadratic approximation in the second variable, which reduces the amount 
of data by two- thirds but increases the evaluation time slightly. 
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i 


2.50E- 


01 


3.50E- 


01 


4.60E- 


•01 


5.60E- 


•01 


7.00E- 


•01 




8.40E- 


•01 


9.80E- 


01 


1.14E 


00 


1.29E 


00 


1.45E 


00 




1.61E 


00 


1.80E 


00 


1.99E 


00 


2. 19E 


00 


2.43E 


00 




2.68E 


00 


2.91E 


00 


3. 14E 


00 


3.36E 


00 


3.56E 


00 




3.77E 


00 


















1 


0.0 




l.OOE- 


01 


2.10E- 


•01 


3.30E- 


01 


4.60E- 


•01 




5.90E- 


01 


7.10E- 


01 


8.40E- 


•01 


1 .OOE 


00 


1.17E 


00 




1.34E 


00 


1.53E 


00 


1.72E 


00 


1.92E 


00 


2.13E 


00 




2. 35E 


00 


2.58E 


00 


2. 83E 


00 


3. 08E 


00 


3.32E 


00 




3.58E 


00 


















1 


-2.50E- 


01 


-1.40E- 


01 


-2.00E- 


•02 


1. lOE- 


•01 


2.40E- 


•01 




3.70E- 


•01 


5.20E- 


•01 


6.70E- 


•01 


8.30E- 


•01 


9.80E- 


•01 




1.15E 


00 


1.33E 


00 


1.50E 


00 


1.70E 


00 


1.89E 


00 




2.10E 


00 


2.31E 


00 


2.54E 


00 


2. 79E 


00 


3.03E 


00 




3.27E 


go 


















1 


-3.70E- 


01 


-2.60E- 


01 


-1 .40E- 


•01 


-1 .OOE- 


•02 


1.20E- 


•01 




2.60E- 


•01 


4.00E- 


01 


5.50E- 


•01 


6. 90E- 


01 


8.60E- 


•01 




I.OIE 


00 


1. 18E 


00 


1.37E 


00 


1.55E 


00 


1.76E 


00 




1.95E 


00 


2. 18E 


00 


2.39E 


00 


2.64E 


00 


2.91E 


00 



3.18E 00 
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TABLE 2 CONTINUED 



1 


-4.&0E-01 


-3.70E-01 


-2.50E- 


■01 


-1.30E- 


■01 


0.0 






1.30E-01 


2.80E-01 


4.20E- 


•01 


5.80E- 


01 


7.40E^ 


-01 




9.00E-01 


1.07E 00 


1.25E 


00 


1.44E 


00 


1.65E 


00 




1.87E 00 
2.95E 00 


2.09E 00 


2.31E 


00 


2. 53E 


00 


2.74E 


00 


1 


-6.20E-01 


-5.00E-01 


-3.80E- 


•01 


-2.50E- 


■01 


•-1 .lOE- 


-01 




2.00E-02 


1.60E-01 


3.10E- 


01 


4.60E- 


■01 


6.00E- 


-01 




7.60E-01 


9.30E-01 


l.llE 


00 


1.30E 


00 


1.50E 


00 




1.71E 00 
2.79E 00 


1.92E 00 


2.15E 


00 


2.35E 


00 


2.58E 


00 


1 


5.70E-01 


6.50E-01 


7.60E- 


■01 


8.70E- 


■01 


9.90E 


-01 




1.13E 00 


1.27E 00 


1.45E 


00 


1.64E 


00 


1.84E 


00 




2.06E 00 


2.27E 00 


2.47E 


00 


2.66E 


00 


2.85E 


00 




3.03E 00 
3.84E 00 


3.21E 00 


3.37E 


00 


3.53E 


00 


3.68E 


00 


1 


4.70E-01 


5.70E-01 


6.90E- 


■01 


8.00E- 


■01 


9.20E 


-01 




1.04E 00 


1.16E 00 


1.30E 


00 


1.47E 


00 


1.65E 


00 




1.86E 00 


2.10E 00 


2.30E 


00 


2.54E 


00 


2.73E 


00 




2.92E 00 
3.84E 00 


3. HE 00 


3.3CE 


00 


3.48E 


00 


3.66E 


00 


1 


3.50E-01 


4.40E-01 


5.50E- 


■01 


6.60E- 


■01 


7.60E- 


-01 




8.90E-01 


1.02E 00 


1 .15E 


00 


1.31E 


00 


1.47E 


00 




1.65E 00 


1.86E 00 


2.06E 


00 


2.27E 


00 


2.48E 


00 




2.71E 00 
3.84E 00 


2.93E 00 


3.14E 


00 


3.37E 


00 


3.61E 


00 


1 


2.10E-01 


3.10E-01 


4.10E- 


■01 


5.20E- 


■01 


6.30E^ 


-01 




7.50E-01 


8.70E-01 


I.OIE 


00 


1.14E 


00 


1.31E 


00 




1.47E 00 


1.66E 00 


1.91E 


00 


2. 13E 


00 


2.35E 


00 




2.57E 00 
3.84E 00 


2.81E 00 


3.05E 


00 


3.29E 


00 


3.56E 


00 


1 


0.0 


8.00E-02 


1.80E- 


■01 


2.80E- 


01 


4.00E- 


-01 




5.10E-01 


6.30E-01 


7.50E- 


■01 


8.90E- 


■01 


1.03E 


00 




1.19E 00 


1.35E 00 


1.53E 


00 


1.73E 


00 


1.94E 


00 




2.17E 00 
3.52E 00 


2.40E 00 


2.65E 


00 


2.S3E 


00 


3.22E 


00 


1 


-1.90E-01 


-9.00E-02 


l.OOE- 


■02 


1.20E- 


•01 


2.30E 


-01 




3.40E-01 


4.60E-01 


5.90E- 


■01 


7.30E-01 


8.70E^ 


-01 




1.02E 00 


1.17E 00 


1.35E 


00 


1. 53E 


00 


1.73E 


00 




1.93E 00 
3.08E 00 


2.15E 00 


2.36E 


00 


2. 58E 


00 


2.83E 


00 


1 


-3.00E-01 


-2.00E-01 


-l.OOE- 


■01 


1 .OOE- 


■02 


1.20E 


-01 




2.40E-01 


3.70E-01 


5.00E- 


01 


6.40E- 


■01 


7.90E^ 


-01 




9.30E-01 


1.09E 00 


1.25E 


00 


1.42E 


00 


1.62E 


00 




1.81E 00 
2.93E 00 


2.03E 00 


2.24E 


00 


2.46E 


00 


2.70E 


00 


1 


-4.10E-01 


-3.00E-01 


-1.90E- 


■01 


-l.OOE- 


■01 


2.00E- 


-02 




1.40E-01 


2.70E-01 


3.90E- 


■01 


5.20E- 


■01 


6.70E 


-01 




8.20E-01 


9.70E-01 


1.14E 


00 


1.22E 


00 


1.51E 


00 




1.70E 00 
2.81E 00 


1.92E 00 


2.13E 


00 


2. 35E 


00 


2.58E 


00 


1 


-5.20E-01 


-4.30E-01 


-3.30E- 


■01 


-2.30E- 


■01 


-I.IOE- 


-01 




l.OOE-02 


1.20E-01 


2.60E- 


■01 


4. COE-01 


5.50E- 


-01 




7.00E-01 


8.60E-01 


1.03E 


00 


1.20E 


00 


1.40E 


00 




1.59E 00 
2.64E 00 


1.78E 00 


1.99E 


00 


2.20E 


00 


2.41E 


00 


1 


5.20E-01 


6.20E-01 


7.00E- 


■01 


8.10E- 


■01 


9.40E- 


-01 




1.08E 00 


1.22E 00 


1.39E 


00 


1.57E 


00 


1.75E 


00 




1.96E 00 


2.18E 00 


2.40E 


00 


2.59E 


00 


2.81E 


00 




3.00E 00 
3.86E 00 


3.19E 00 


3.37E 


00 


3.56E 


00 


3.71E 


00 


1 


3.90E-01 


5.00E-01 


6.10E- 


■01 


7.30E- 


■01 


8.30E 


-01 




9.60E-01 


l.lOE 00 


1.24E 


00 


1 .40E 


00 


1.57E 


00 




1.76E 00 


1.97E 00 


2.19E 


00 


2.42E 


00 


2.64E 


00 




2.86E 00 
3.86E 00 


3.06E 00 


3.27E 


00 


3.47E 


00 


3.67E 


00 


1 


2.70E-01 


3.60E-01 


4.60E- 


■01 


5 .60E- 


■01 


6.80E- 


-01 




8.10E-01 


9.40E-01 


1.09E 


00 


1 .24E 


00 


1.40E 


00 




1.58E 00 


1.78E 00 


1.97E 


00 


2. 19E 


00 


2.41E 


00 




2.65E 00 


2.89E 00 


3.13E 


00 


3.36E 


00 


3.59E 


00 



3.80E 00 
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^811 

5 8 11 

6 8 11 

7 8 11 

8 8 11 

9 8 11 



1.90E-01 
7.20E-01 
1.45E 00 
2.48E 00 
3.80E 00 
0.0 

5.00E-01 
1.20E 00 
2.16E 00 
3.53E 00 
-1.80E-01 
3.60E-01 
1.03E 00 
1.96E 00 
3.18E 00 
-2.70E-01 
2.60E-01 
9.30E-01 
1.87E 00 
3.00E 00 
-3.80E-01 
1.40E-01 
8.40E-01 
1.76E 00 
2.89E 00 
-5.00E-01 
4.00E-02 
7.30E-01 
1.59E 00 
2.70E 00 



2.90E-01 
8.50E-01 
1.62E 00 
2.71E 00 

8.00E-02 
6.30E-01 
1.37E 00 
2.39E 00 

-9.00E-02 
4.90E-01 
1.18E 00 
2.18E 00 

-1.60E-01 
4.00E-01 
1.08E 00 
2.08E 00 

-2.80E-01 
2.80E-01 
l.OOE 00 
1.95E 00 

-4.20E-01 
1.80E-01 
8.70E-G1 
1.79E 00 



3.80E-01 
9.80E-01 
1.82E 00 
2.95E 00 

1.70E-01 
7.70E-01 
1.54E 00 
2.65E 00 

2.00E-02 
6.30E-01 
1.36E 00 
2.42E 00 

-6.00E-02 
5.30E-01 
1.26E 00 
2.30E 00 

-1.90E-01 
4.30E-01 
1.17E 00 
2.18E 00 

-3.10E-01 
3.10E-01 
1.03E 00 
2.01E 00 



4. SOE-01 

1. 12E 00 

2. C3E 00 
3.21E 00 

2.70E-01 
9. lOE-01 
1.72E 00 
2.92E 00 

1.20E-01 
7. 70E-01 
1.56E 00 
2.66E 00 

5. OOE-02 
6.60E-01 
1.45E 00 
2.53E 00 

-9.00E-02 
5.70E-01 
1.36E 00 
2.40E 00 

-2.10E-01 
4.60E-01 
1.20E 00 
2.23E 00 



6.10E-01 
1.28E 00 
2.25E 00 
3.49E 00 

3.90E-01 
1.05E 00 
1.92E 00 
3.21E 00 

2.40E-01 
9.00E-01 
1.76E 00 
2.91E 00 

1.50E-01 
7.90E-01 
1.66E 00 
2.77E 00 

3.00E-02 
7 .lOE-01 
1.55E OC 
2.65E OC 

-9.00E-02 
6.00E-01 
1.38E 00 
2.47E 00 
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THE 21 VALUES OF THE ITH VARIABLE FELLOW 











0.0 




2.00E 


00 


4.00E 


00 


6.00E 


00 


8.00E 


00 










l.OOE 


01 


1.20E 


01 


1.40E 


01 


1.60E 


01 


1.80E 


01 










2.00E 


01 


2.20E 


01 


2.40E 


01 


2.60E 


01 


2.80E 


01 










3.00E 


01 


3.20E 


01 


3.40E 


01 


3.60E 


01 


3.80E 


01 










4.00E 


01 
























THE 9 VALUES OF THE 2TH VARIABLE FELLOW 














2.00E 


01 


1.60E 


01 


1.20E 


01 


8.00E 


00 


O 

• 

o 












-8.00E 


00 


-1.20E 


01 


-1.60E 


01 


-2.00E 


01 












THE 8 VALUES OF THE 3TH VARIABLE FELLOW 














1.30E 


00 


1.60E 


00 


2.00E 


00 


2.30E 


00 


2.60E 


00 










3.00E 


00 


4.00E 


00 


4.60E 


00 










12 


13 


14 


15 


11 = 1 




11=2 




ETC. 












1 


1 


1 


1 


6.40E-01 


7.90E- 


01 


9.60E- 


■01 


l.lOE 


00 


1.29E 


00 










1.46E 


00 


1.64E 


00 


1.84E 


00 


2.03E 


00 


2.22E 


00 










2.42E 


00 


2.63E 


00 


2.83E 


00 


3.C5E 


00 


3.28E 


00 










3.49E 


00 


3.72E 


00 


3.94E 


00 


4.18E 


00 


4.42E 


00 










4.64E 


00 


















2 


1 


1 


1 


5.70E- 


01 


7.40E- 


01 


9.10E- 


•01 


l.lOE 


00 


1.29E 


00 










1.47E 


00 


1.63E 


00 


1.83E 


00 


2.01E 


00 


2.22E 


00 










2.42E 


00 


2.63E 


00 


2.83E 


00 


3.C5E 


00 


3.28E 


00 










3.50E 


00 


3.73E 


00 


3.95E 


00 


4.18E 


00 


4.41E 


00 










4.64E 


00 


















3 


1 


1 


1 


4.70E- 


■01 


6.60E- 


■01 


8.60E- 


•01 


1.05E 


00 


1.24E 


00 










1.44E 


00 


1.63E 


00 


1.83E 


00 


2.02E 


00 


2.22E 


00 










2.42E 


00 


2.63E 


00 


2.85E 


00 


3.C5E 


00 


3.27E 


00 










3.50E 


00 


3.72E 


00 


3.95E 


00 


4.19E 


00 


4.43E 


00 










4.64E 


00 


















4 


1 


1 


1 


3.10E- 


•01 


5.20E- 


•01 


7.20E- 


•01 


9.40E- 


-01 


1.15E 


00 










1.36E 


00 


1.56E 


00 


1.76E 


00 


1.95E 


00 


2.15E 


00 










2.35E 


00 


2.55E 


00 


2.76E 


00 


2.S6E 


00 


3.15E 


00 










3.36E 


00 


3.58E 


00 


3.82E 


00 


4.05E 


00 


4.31E 


00 










4.57E 


00 


















5 


1 


1 


1 


0.0 




2.10E- 


01 


4.50E- 


•01 


6.80E-01 


9.00E- 


•01 










1.12E 


00 


1.35E 


00 


1.58E 


00 


1.79E 


00 


2.01E 


00 










2.23E 


00 


2.46E 


00 


2.69E 


00 


2.91E 


00 


3.13E 


00 










3.36E 


00 


3.59E 


00 


3.81E 


00 


4.05E 


00 


4.27E 


00 










4.48E 


00 


















6 


1 


1 


1 


-3.10E- 


•01 


-l.OOE- 


01 


I.IOE- 


•01 


3.30E- 


■01 


5.60E- 


-01 










7.80E- 


•01 


l.OOE 


00 


1.23E 


00 


1.46E 


00 


1.69E 


00 










1.92E 


00 


2.15E 


00 


2.39E 


00 


2.63E 


00 


2.88E 


00 










3.13E 


00 


3.36E 


00 


3.61E 


00 


3.85E 


00 


4.09E 


00 










4.34E 


00 


















7 


1 


1 


1 


-4.60E- 


•01 


-2.70E- 


01 


-7.00E- 


•02 


1.30E- 


■01 


3.50E- 


■01 










5.70E- 


•01 


7.90E- 


01 


l.OOE 


00 


1.25E 


00 


1.48E 


00 










1.71E 


00 


1.96E 


00 


2.19E 


00 


2.43E 


00 


2.68E 


00 










2.92E 


00 


3.16E 


00 


3.42E 


00 


3.66E 


00 


3.92E 


00 










4.18E 


00 


















8 


1 


1 


1 


-5.50E- 


•01 


-3.80E- 


01 


-1.80E- 


•01 


l.OOE- 


■02 


2.00E- 


■01 










4.10E- 


•01 


6.20E- 


01 


8 .50E- 


•01 


1.C7E 


00 


1.32E 


00 










1.55E 


00 


1.80E 


00 


2.05E 


00 


2.29E 


00 


2.54E 


00 










2.79E 


00 


3.03E 


00 


3.26E 


00 


3.54E 


00 


3.78E 


00 










4.05E 


00 


















9 


1 


1 


1 


-6.50E- 


•01 


-4.50E- 


•01 


-2.60E- 


•01 


-7.00E- 


■02 


1.20E- 


■01 










3.40E- 


•01 


5.50E- 


01 


7.70E'01 


l.COE 


00 


1.22E 


00 










1.46E 


00 


1.69E 


00 


1.S3E 


00 


2.17E 


00 


2.41E 


00 










2.66E 


00 


2.90E 


00 


3.16E 


00 


3.40E 


00 


3.65E 


00 










3.91E 


00 


















1 


2 


1 


1 


5.40E- 


•01 


7.30E- 


•01 


9.30E- 


■01 


1.12E 


00 


1.28E 


00 










1.45E 


00 


1.62E 


00 


1.78E 


00 


1.94E 


00 


2.12E 


00 










2.30E 


00 


2.49E 


00 


2.69E 


00 


2.91E 


00 


3.14E 


00 










3.37E 


00 


3.62E 


00 


3.86E 


00 


4. 12E 


00 


4.37E 


00 










4.65E 


00 
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6 

7 

8 

9 

1 

2 

3 
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5 
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2 



2 



2 



2 



2 



2 



3 



3 



3 



3 



3 
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1 


1 


5.10E- 


01 


6.80E- 


01 


8.50E- 


•01 


1 .C3E 


00 


1.21E 


00 






1.38E 


00 


1.55E 


00 


1.73E 


00 


1 .91E 


00 


2.08E 


00 






2.28E 


00 


2.48E 


00 


2.68E 


00 


2.91E 


00 


3.13E 


00 






3.37E 


00 


3.61E 


00 


3.87E 


00 


4. 13E 


00 


4.38E 


00 






4.65E 


00 


















1 


1 


3.90E- 


■01 


5.70E- 


01 


7.50E- 


•01 


9.10E- 


•01 


1.09E 


00 






1.28E 


00 


1.46E 


00 


1 .67E 


00 


1.85E 


00 


2.04E 


00 






2.25E 


00 


2.47E 


00 


2.69E 


00 


2.91E 


00 


3.14E 


00 






3.37E 


00 


3.62E 


00 


3.87E 


00 


4.14E 


00 


4.38E 


00 






4.65E 


00 


















1 


1 


2.60E- 


•01 


4.40E- 


01 


6 .60E- 


■01 


6.40E- 


•01 


1.04E 


00 






1.22E 


00 


1.42E 


00 


1 .60E 


00 


1.77E 


00 


1.95E 


00 






2.16E 


00 


2.38E 


00 


2.60E 


00 


2.83E 


00 


3 .06E 


00 






3.30E 
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-1, 


,47E 


00 


-1.54E 


00 


-1 


. 60E 


00 


“1 • 66E 


00 


-1 


.73E 


00 










-1, 


.86E 


00 
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TABLE 5 



THE 21 VALUES OF THE ITH VARIABLE FOLLOW 











0.0 




2.00E 


00 


4. COE 


00 


6. COE 


00 


8. OOE 


00 










l.OOE 


01 


1.20E 


01 


1.40E 


01 


1 .60E 


01 


1.80E 


01 










2.00E 


01 


2.20E 


01 


2.40E 


01 


2.60E 


01 


2.80E 


01 










3.00E 


01 


3.20E 


01 


3.40E 


01 


3.60E 


01 


3.80E 


01 










4.00E 


01 
























THE 9 VALUES OF THE 2TH VARIABLE FOLLOW 














2.00E 


01 


1.60E 


01 


1.20E 


01 


8 .OOE 


00 


0.0 












-8.00E 


00 


-1.20E 


01 


-1.60E 


01 


-2.00E 


01 












THE 8 VALUES OF THE 3TH VARIABLE FOLLOW 














1.30E 


00 


1.60E 


00 


2.00E 


00 


2.30E 


00 


2.60E 


00 










3.00E 


00 


4.00E 


00 


4.60E 


00 










12 


13 


14 


15 


11 = 1 




11=2 




ETC. 












1 


1 


1 


1 


3.06E 


00 


3.13E 


00 


3.09E 


00 


2.98E 


00 


2.76E 


00 










2.26E 


00 


1.69E 


00 


1.23E 


00 


8.50E- 


•01 


5.00E- 


-01 










1.40E- 


'01 


-2.30E-01 


-5.90E- 


•01 


-9.50E- 


-01 


-1.27E 


00 










-1.50E 


00 


-1.72E 


00 


-1.93E 


00 


-2. HE 


00 


-2.24E 


00 










-2.35E 


00 


















2 


1 


1 


1 


2.61E 


00 


2.64E 


00 


2.57E 


00 


2.38E 


00 


2.06E 


00 










1.57E 


00 


l.lOE 


00 


7.30E- 


•01 


4.20E- 


•01 


l.OOE- 


-01 










-2.40E- 


•01 


-5.70E- 


•01 


-8.80E-01 


-1.22E 


00 


-1.49E 


00 










-1.72E 


00 


-1.93E 


00 


-2.10E 


00 


-2.27E 


00 


-2.41E 


00 










-2.51E 


00 


















3 


1 


1 


1 


1.98E 


00 


1.94E 


00 


1.82E 


00 


1.60E 


00 


1.31E 


00 










9.50E- 


•01 


5.70E- 


01 


2.90E- 


•01 


2.00E- 


-02 


-2.80E- 


-01 










-5.60E- 


•01 


-8.40E- 


•01 


-1.12E 


00 


-1.39E 


00 


-1.62E 


00 










-1.84E 


00 


-2.04E 


00 


-2.21E 


00 


-2.35E 


00 


-2.49E 


00 










-2.59E 


00 


















4 


1 


1 


1 


1.29E 


00 


1.24E 


00 


l.lOE 


00 


9. OOE- 


•01 


6.30E- 


-01 










3.60E- 


01 


I.IOE- 


•01 


-1.40E- 


■01 


-3 .90E-01 


-6.30E- 


-01 










-8.90E- 


•01 


-1.14E 


00 


-1.39E 


00 


-1.63E 


00 


-1.84E 


00 










-2.03E 


00 


-2.22E 


00 


-2.38E 


00 


-2.51E 


00 


-2.63E 


00 










-2.73E 


00 


















5 


1 


1 


1 


0.0 




-1.50E- 


•01 


-2.80E- 


•01 


-4.10E- 


•01 


-5.60E- 


-01 










-7.20E- 


■01 


-8.70E-01 


-1.03E 


00 


-1.22E 


00 


-1.39E 


00 










-1.55E 


00 


-1.73E 


00 


-1.90E 


00 


-2.05E 


00 


-2.20E 


00 










-2.33E 


00 


-2.45E 


oc* 


-2.56E 


00 


-2.64E 


00 


-2.73E 


00 










-2.84E 


00 


















6 


1 


1 


1 


-1.32E 


00 


-1.33E 


00 


-1.36E 


00 


-1.42E 


00 


-1.52E 


00 










-1.68E 


00 


-1.80E 


00 


-1.91E 


00 


-2.01E 


00 


-2. HE 


00 










-2.18E 


00 


-2.27E 


00 


-2.34E 


00 


-2.41E 


00 


-2.47E 


00 










-2.54E 


00 


-2.58E 


00 


-2.65E 


00 


-2.71E 


00 


-2.78E 


00 










-2.85E 


00 


















7 


1 


1 


1 


-1.99E 


00 


-1.93E 


00 


-1.92E 


00 


-1.95E 


00 


'2.03E 


00 










-2.15E 


00 


-2.27E 


00 


-2.36E 


00 


-2.43E 


00 


-2.49E 


00 










-2.53E 


00 


-2.56E 


00 


^2.58E 


00 


-2.62E 


00 


-2.65E 


00 










-2.68E 


00 


-2.73E 


00 


-2.77E 


00 


-2.81E 


00 


-2.87E 


00 










-2.92E 


00 


















8 


1 


1 


1 


-2.63E 


00 


-2.49E 


00 


-2.42E 


00 


-2.42E 


00 


■~2 .47E 


00 










-2.57E 


00 


-2.66E 


00 


-2.74E 


00 


-2.79E 


00 


-2.84E 


00 










-2.87E 


00 


-2.87E 


00 


-2.86E 


00 


-2.85E 


00 


-2.85E 


00 










-2.85E 


00 


-2.87E 


00 


-2.89E 


00 


-2. 91E 


00 


-2.94E 


00 










-2.97E 


00 


















9 


1 


1 


1 


-3.08E 


00 


-2.92E 


00 


-2.83E 


00 


-2.81E 


00 


-2.83E 


00 










-2.92E 


00 


-2.99E 


00 


-3.09E 


00 


-3. 18E 


00 


-3.24E 


00 










-3.28E 


00 


-3.29E 


00 


-3.26E 


00 


-3. 22E 


00 


-3.20E 


00 










-3.16E 


00 


-3.15E 


00 


-3.14E 


00 


-3. 14E 


00 


-3.15E 


00 










-3.15E 


00 


















1 


2 


1 


1 


2.08E 


00 


2.04E 


00 


1.95E 


00 


1.79E 


00 


1.57E 


00 










1.23E 


00 


8.70E- 


•01 


5.40E- 


•01 


2.50E- 


•01 


4.00E- 


-02 










-2.00E- 


01 


-4.30E- 


•01 


-6.50E-01 


-8.80E- 


•01 


-l.llE 


00 










-1.31E 


00 


-1.48E 


00 


-1.67E 


00 


-1.81E 


00 


-1.96E 


00 










-2.08E 


00 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



TABLE 5 CONTINUED 



1 


1, 


,59E 


00 


1 


.52E 


00 


1 


.42E 


00 


1 


.25E 


00 


1 


.05E 


00 




8. 


. OOE- 


-01 


5 


.30E- 


-01 


3 


.OOE- 


-01 


6 


. COE- 


-02 


-1 


.70E- 


-01 




-3< 


,80E- 


-01 


-5 


.60E- 


-01 


-7 


.70E-01 


-9 


.60E- 


-01 


-1 


.17E 


00 




-1, 


,37E 


00 


-1 


. 55E 


00 


-1 


. 72E 


00 


-1 


.87E 


00 


-2 


.02E 


00 




-2. 


.13E 


00 


























1 


1, 


.13E 


00 


1 


.04E 


00 


9 


.20E- 


-01 


7 


. 70E- 


•01 


5 


.70E- 


-01 




3. 


,60E- 


-01 


1 


.40E- 


-01 


-4 


.OOE- 


-02 


-2 


.30E- 


-01 


-4 


.30E- 


-01 




-6. 


.lOE- 


-01 


-8 


.OOE- 


-01 


-9 


.80E- 


-01 


-1 


. 16E 


00 


-1 


.32E 


00 




-1. 


. 50E 


00 


-1 


• 66 E 


00 


-1 


.80E 


00 


-1 


• 95E 


00 


-2 


.06E 


00 




-2. 


.19E 


00 


























1 


7. 


,70E« 


-01 


6 


.50E- 


-01 


5 


.20E- 


-01 


3 


.80E- 


-01 


2 


.30E- 


-01 




7. 


,00E- 


-02 


-1 


. OOE- 


-01 


-2 


.70E-01 


-4 


. 60E- 


-01 


-6 


.40E-01 




-8. 


• 20E- 


-01 


-9 


.90E- 


-01 


-1 


.17E 


00 


-1 


. 34E 


00 


-1 


.50E 


00 




-1, 


» 66E 


00 


-1 


.81E 


00 


-1 


• 96E 


00 


-2 


. C8E 


oc 


-2 


.18E 


00 




-2, 


.29E 


00 


























1 


0< 


>0 




-1 


.40E- 


■01 


-2 


. 70E- 


-01 


-4 


. OOE- 


-01 


-5 


.20E- 


-01 




-6. 


• 50E- 


-01 


-7 


.60E- 


-01 


-8 


. 50E- 


-01 


-9 


. 80E- 


-01 


-1 


.lOE 


00 




-1. 


,21E 


00 


-1 


.32E 


00 


-1 


.45E 


00 


-1 


.58E 


00 


-1 


.72E 


00 




-1. 


,86E 


00 


-1 


.99E 


00 


-2 


. 13E 


00 


-2 


.26E 


00 


-2 


.35E 


00 




-2. 


,44E 


00 


























1 


-7, 


.80E- 


-01 


-8 


.40E- 


-01 


-9 


.20E- 


-01 


-9 


. 90E- 


-01 


-1 


.09E 


00 




-1. 


. 19E 


00 


-1 


.28E 


00 


-1 


.39E 


00 


-1 


. A8E 


00 


-1 


.60E 


00 




-1. 


,70E 


oc 


-1 


. 79E 


00 


-1 


.91E 


00 


-2 


.OlE 


00 


-2 


.12E 


00 




-2, 


,21E 


00 


-2 


.31E 


00 


-2 


.37E 


00 


-2 


. 46E 


00 


-2 


.55E 


oc 




-2, 


> 6 1 E 


00 


























1 


-1, 


,13E 


00 


-1 


.22E 


00 


-1 


.29E 


00 


-1 


.36E 


00 


-1 


.42E 


00 




-1. 


,51E 


00 


-1 


.59E 


00 


-1 


.73E 


00 


-1 


. 84E 


00 


-1 


.94E 


00 




-2. 


,03E 


00 


-2 


.12E 


00 


-2 


. 18E 


00 


-2 


.27E 


00 


-2 


.33E 


00 




-2. 


,42E 


00 


-2 


.48E 


00 


-2 


.56E 


00 


-2 


• 64E 


00 


-2 


.71E 


00 




-2. 


► 79E 


00 


























1 


-1. 


,60E 


00 


-1 


.62E 


00 


-1 


.62E 


00 


-1 


• 66 E 


00 


-1 


.71E 


00 




-1< 


,77E 


00 


-1 


.85E 


00 


-1 


.97E 


00 


-2 


.C9E 


00 


-2 


.17E 


00 




-2, 


,27E 


00 


-2 


.35E 


00 


-2 


.41E 


00 


-2 


*49E 


00 


-2 


.55E 


00 




-2, 


.62E 


00 


-2 


.69E 


00 


-2 


.75E 


00 


-2 


.81E 


00 


-2 


• 86E 


00 




'2. 


► 91E 


00 


























1 


-2. 


,08E 


00 


-2 


.06E 


00 


-2 


.04E 


00 


-1 


.98E 


00 


-1 


.98E 


00 




-2. 


.02E 


00 


-2 


.lOE 


00 


-2 


.21E 


00 


-2 


.35E 


00 


-2 


.46E 


00 




-2. 


► 54E 


00 


-2 


.61E 


00 


-2 


.68E 


00 


-2 


.73E 


00 


-2 


.78E 


00 




-2, 


,83E 


00 


-2 


.88E 


00 


-2 


.94E 


00 


-2 


.98E 


00 


-3 


.04E 


00 




-3< 


► 09E 


00 


























1 


1. 


,16E 


00 


1 


. 13E 


00 


1 


.06E 


00 


9 


. 50E- 


-01 


8 


.OOE- 


-01 




6. 


.30E-01 


4 


.lOE- 


-01 


2 


.50E- 


-01 


1 


.OOE- 


-01 


-4 


.OOE- 


-02 






,70E- 


-01 


-3 


.50E- 


-01 


-5 


.lOE- 


-D1 


-7 


.lOE- 


-01 


-8 


.80E- 


-01 




-1. 


, 09E 


00 


-1 


.33E 


00 


-1 


.53E 


00 


-1 


.78E 


00 


-2 


.03E 


00 




-2. 


,29E 


00 


























1 


9. 


.30E-01 


8 


.70E- 


-01 


7 


.90E- 


-01 


6 


.70E- 


■01 


5 


.40E- 


-01 




3. 


.80E- 


-01 


2 


.lOE- 


-01 


9 


.OOE- 


-02 


-4 


. OOE- 


-02 


-1 


.60E- 


-01 




-2. 


,80E- 


-01 


-4 


.40E- 


-01 


-6 


.lOE- 


-01 


-8 


. lOE- 


-01 


“9 


.80E- 


-01 




-1, 


,21E 


00 


-1 


.42E 


00 


-1 


.63E 


00 


-1 


.89E 


00 


-2 


.14E 


00 




-2, 


.36E 


00 


























1 


6. 


,80E- 


-01 


6 


.lOE- 


-01 


5 


.lOE-Ol 


3 


.90E- 


-01 


2 


.70E- 


-01 




1. 


.50E-01 


2 


.OOE- 


-02 


-8 


.OOE-02 


-1 


. 50E- 


-01 


-2 


.40E- 


-01 




-3, 


,70E- 


-01 


-5 


.lOE- 


-01 


-6 


.90E- 


-01 


-8 


.90E-01 


-1 


.09E 


00 




-1, 


,31E 


00 


-1 


.52E 


00 


-1 


.76E 


00 


-2 


.OOE 


00 


-2 


.25E 


oc 




-2. 


,50E 


00 


























1 


4, 


,90E> 


-01 


3 


.80E- 


-01 


2 


.40E- 


-01 


1 


. lOE-Ol 


-2 


.OOE- 


-02 




-1. 


.lOE-01 


-1 


.50E- 


-01 


-2 


.30E-01 


-3 


. lOE- 


-01 


-4 


.20E- 


-01 




-5. 


. 50E- 


-01 


-7 


. OOE- 


-01 


-8 


.60E- 


-01 


-1 


.C6E 


00 


-1 


.26E 


00 




-1. 


,48E 


00 


-1 


.69E 


00 


-1 


.93E 


00 


-2 


. 15E 


00 


-2 


.40E 


00 




-2. 


,63E 


00 


























1 


0. 


.0 




-1 


.20E- 


-01 


-2 


. 50E- 


-01 


-3 


. 50E- 


-01 


-4 


.40E- 


-01 




-5. 


,40E- 


-01 


-6 


.iOE- 


-01 


-6 


.70E- 


-01 


-7 


. 50E- 


-01 


-*•8 


.40E- 


-01 




-9, 


,30E- 


-01 


-1 


.05E 


00 


-1 


.17E 


00 


-1 


.35E 


00 


-1 


.55E 


00 




-1, 


.77E 


00 


-1 


.98E 


00 


-2 


.21E 


00 


-2 


. A5E 


00 


-2 


.68E 


00 




-2. 


► 86E 


00 


























1 


-4, 


,70E- 


-01 


-5 


.50E- 


-01 


-6 


.30E-01 


-7 


. lOE- 


-01 


-8 


.OOE- 


-01 




-8, 


.90E- 


-01 


-9 


.80E- 


-01 


-1 


.06E 


00 


-1 


. 13E 


00 


-1 


.21E 


00 




-1, 


,31E 


00 


-1 


.43E 


00 


-1 


. 58E 


00 


-1 


. 77E 


00 


-1 


.98E 


00 




-2, 


► 21E 


00 


-2 


.40E 


00 


-2 


.61E 


00 


-2 


.82E 


00 


-3 


.OOE 


00 




-3. 


,16E 


00 
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-1 . lOE-01 


-1.60E-01 






-2.30E-01 


-3.10E-01 


-4.00E-01 


-5.C0E-01 


-6.00E-01 






-7.10E-01 


-8.60E-«01 


-1.03E 00 


-1.22E 00 


-1.43E 00 






-1.65E 00 
-2.99E 00 


-1.91E 00 


-2.17E 00 


-2.42E 00 


-2.71E 00 


1 


1 


-3.00E-02 


-6.00E-02 


-l.lOE-01 


-1.60E-01 


-2.10E-01 






-2.90E-01 


-3.80E-01 


-4.90E-01 


-6.COE-01 


-7.30E-01 






-8.70E-0I 


-1.02E 00 


-1.21E OO 


-1.41E 00 


-1.65E 00 






-1.89E 00 
-3. HE 00 


-2.14E 00 


-2.39E 00 


-2.62E 00 


-2.87E 00 


1 


1 


-7.00E-02 


-9.00E-02 


-1.30E-01 


-1.80E-01 


-2.40E-01 






-3.20E-01 


-4.20E-01 


-5.50E-01 


-6.80E-01 


-8.20E-01 






-9.70E-01 


-1.15E 00 


-1.34E 00 


-1.58E 00 


-1.82E 00 






-2.10E 00 
-3.27E 00 


-2.35E 00 


-2.60E 00 


-2.84E 00 


-3.06E 00 


1 


1 


-9.00E-02 


-l.lOE-01 


-1 .60E-01 


-2. 20E-01 


-2.90E-01 






-3.80E-01 


-4.80E-01 


-6.20E-01 


-7.70E-01 


-9.20E-01 






-1.08E 00 


-1.26E 00 


-1.49E 00 


-1.76E 00 


-2.05E 00 






-2.37E 00 
-3.50E 00 


“2.66E 00 


-2.92E 00 


-3.16E 00 


-3.35E 00 


1 


1 


8.00E-02 


l.lOE-01 


1.40E-01 


1.70E-01 


2.00E-01 






2.30E-01 


2.60E-01 


2.80E-01 


2.90E-01 


3 .OOE-01 






2.90E-01 


2.70E-01 


2.70E-01 


2.50E-01 


2.10E-01 






1.70E-01 

-8.60E-01 


7.00E-02 


-7.00E-02 


-2.80E-01 


-5.40E-01 


1 


1 


6.00E-02 


7.00E-02 


9.00E-02 


l.lOE-01 


1.30E-01 






1.60E-01 


1.80E-01 


2.00E-01 


2.10E-01 


2.10E-01 






2.10E-01 


2.006-01 


1.80E-01 


1.60E-01 


1.30E-01 






8.00E-02 
-1.08E 00 


0.0 


-1.60E-01 


-3.90E-01 


-7.00E-01 


1 


1 


3.00E-02 


3.00E-02 


5.00E-02 


5.C0E-02 


6.00E-02 






7.00E-02 


9.00E-02 


l.lOE-01 


1 . lOE-01 


l.lOE-01 






l.OOE-01 


9.00E-02 


9.00E-02 


8. OOE-02 


6 .OOE— 02 






2.00E-02 
-1.28E 00 


-7.00E-02 


-2.50E-01 


-5. lOE-01 


-8.50E-01 



[34] 



TABLE 5 CONTINUED 



^811 

5 8 11 

6 8 11 

7 8 11 

8 8 11 

9 8 11 



2.00E-02 
l.OOE-02 
'-2.00E-02 
-2.40E-01 
-1.46E 00 
2.00E-02 
-5.00E-02 
-3.80E-01 
-8.70E-01 
-2.14E 00 
0.0 

-1.50E-01 
-7.30E-01 
-1.65E 00 
-2.78E 00 
-2.00E-02 
-1.80E-01 
-8. lOE-01 
-1.93E 00 
-3.26E 90 
-l.OOE-02 
-2.20E-01 
-9.60E-01 
-2.20E 00 
-3.40E 00 
-4.00E-02 
-2.60E-01 
-1.02E 00 
-2.37E 00 
-3.50E 00 



2.00E-02 

l.OOE-02 

-2.00E-02 

-3.80E-01 

2.00E-02 
-l.OOE-01 
-4.60E-01 
-I.OIE 00 

0.0 

-2.30E-01 
-8.80E-01 
-1.88E 00 

-2.00E-02 
-2.60E-01 
-l.OOE 00 
-2.15E 00 

-3.00E-02 
-3.20E-01 
-1.19E 00 
-2.44E 00 

-4.00E-02 
-3.50E-01 
-1.25E 00 
-2.65E 00 



2.C0E-02 

l.OOE-02 

-2.00E-02 

-5.90E-01 

2.00E-02 
-1.60E-01 
-5.50E-01 
-1.20E 00 

-1.00E'02 
-3.30E-01 
-1.05E 00 
-2.10E 00 

-4.00E-02 
-3.70E-01 
-1.21E 00 
-2.41E 00 

-5.00E-02 
-4.10E-01 
-1.42E 00 
-2.68E 00 

-7.00E-02 
-4.70E-01 
-1.52E 00 
-2.90E 00 



2.C0E-02 
1 .COE-02 
-4.C0E-02 
-8.30E-01 

l.COE-02 
-2.20E-01 
-6.A0E-01 
-1.A5E 00 

-3.00E-02 
-4. 50E-01 
-1.23E 00 



-7.C0E-02 
-5.C0E-01 
-1.A4E 00 
-2.68E 00 

-8.00E-02 
-5.S0E-01 
-1.68E 00 
-2.94E 00 

-l.lOE-01 
-6.40E-01 
-1.7SE 00 
-3.12E 00 



2.00E-02 

0.0 

-1 .lOE-01 
-1.12E 00 

-1.0' E-02 
-2.,9CE-01 
-7.40E-01 
-1.75E 00 

-9.00E-02 
-5.90E-01 
-1.43E 00 
-2.57E 00 

-1.20E-01 
-6.50E-01 
-1.68E 00 
-2.96E 00 

-1.50E-01 
-7 .70E-01 
-1.93E 00 
-3.19E 00 

-1.80E-01 
-8.20E-01 
-2.08E 00 
-3.33E 00 
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TABLE 6 



THE 13 VALUES OF TH 

0.0 2.00E 

l.OOE 01 1.20E 

2.00E 01 2.20E 

THE 9 VALUES OF TH 

-2.00E 01 -1.60E 
8.00E 00 1.20E 

THE 10 VALUES OF TH 

8.00E-01 l.OOE 
2.30E 00 2.60E 



2 


13 


14 


15 


11 = 1 




11 = 2 


1 


1 


1 


1 


-2.37E 


00 


-2.26E 










-2.58E 


00 


-2.70E 










-2.97E 


00 


-2.98E 


2 


1 


1 


1 


-2.12E 


00 


-2. 08E 










-2.31E 


00 


-2.47E 










-2.80E 


00 


-2.81E 


3 


1 


1 


1 


-2.00E 


00 


-1.95E 










-2.22E 


00 


-2.38E 










-2.70E 


00 


-2.73E 


4 


1 


1 


1 


-1.30E 


00 


-1.20E 










-1.63E 


00 


-1.89E 










-2.48E 


00 


-2.53E 


5 


1 


1 


1 


0.0 




—2 . OOE— 










-7.00E- 


-01 


-1.08E 










-2.29E 


00 


-2.41E 


6 


1 


1 


1 


1.30E 


00 


1.32E 










-3.10E-01 


-7.00E- 










-2.05E 


00 


-2.25E 


7 


1 


1 


1 


1.99E 


00 


1.99E 










1.20E-01 


-3.O0E- 










-1.86E 


00 


-2.12E 


8 


1 


1 


1 


2.19E 


00 


2.20E 










2.90E- 


-01 


-2.00E- 










-1.74E 


00 


-2.02E 


9 


1 


1 


1 


2.40E 


00 


2.38E 










4.50E- 


-01 


0.0 










-1.646 


00 


-1.94E 


1 


2 


1 


1 


-2.98E 


00 


-2.86E 










-2.96E 


00 


-3.12E 










-3.77E 


00 


-3.89E 


2 


2 


1 


1 


-2.44E 


00 


-2.39E 










-2.72E 


00 


-2.92E 










-3.56E 


00 


-3.69E 


3 


2 


1 


1 


-2.06E 


00 


-1.90E 










-2. 16E 


00 


-2.32E 










-3.01E 


00 


-3.19E 


4 


2 


1 


1 


-1.30E 


00 


-1.23E 










-1.91E 


00 


-2. HE 










-2.89E 


00 


-3.08E 


5 


2 


1 


1 


0.0 




-7.00E- 










-8.00E- 


-01 


-1.12E 










-2.27E 


00 


-2.50E 


6 


2 


1 


1 


1.34E 


00 


1.37E 










-5.00E-02 


-5.00E- 










-1.96E 


00 


-2.22E 


7 


2 


1 


1 


2.06E 


00 


2.04E 










3.00E- 


-01 


-2.00E- 










-1.75E 


00 


-2.05E 


8 


2 


1 


1 


2.41E 


00 


2.39E 










4.70E- 


-01 


-2.00E- 










-1.64E 


00 


-1.96E 


9 


2 


1 


1 


3.00E 


00 


2.98E 










8.50E- 


-01 


2.50E- 










-1.43E 


00 


-1.74E 



I36J 



ITH VARIABLE FOLLOW 



4. OOE 


00 


6. COE 


00 


8. COE 


00 


1.40E 


01 


1.60E 


01 


1.80E 


01 


2.40E 


01 










2TH VARIABLE FCLLOW 






-1.20E 


01 


-8. OOE 


00 


o 

• 

o 




1.60E 


01 


2. OOE 


01 






3TH VARIABLE FCLLOW 






1 .30E 


00 


1.60E 


00 


2. OOE 


00 


3. OOE 


00 


4. OOE 


00 


4.60E 


00 


ETC. 












-2.27E 


00 


-2.32E 


00 


-2.44E 


00 


-2.81E 


00 


-2.89E 


00 


-2.94E 


00 


-2.94E 


00 










-2.04E 


00 


-2.08E 


00 


-2.16E 


00 


-2.60E 


00 


-2. TOE 


00 


-2.77E 


00 


-2.81E 


00 










-1.94E 


00 


-1.S7E 


00 


-2.07E 


00 


-2.50E 


00 


-2.60E 


00 


-2.67E 


00 


-2.74E 


00 










-1.19E 


00 


-1.25E 


00 


-1.41E 


00 


-2.09E 


00 


-2.24E 


00 


-2.39E 


00 


-2.57E 


00 










-1 . lOE- 


-01 


-2.5QE- 


-01 


-4.20E- 


•01 


-1.47E 


00 


-1.80E 


00 


-2.07E 


00 


-2.45E 


00 










1.09E 


00 


5.C0E- 


•01 


4.00E- 


■02 


-1.09E 


00 


-1.41E 


00 


-1.72E 


00 


-2.34E 


00 










1.73E 


00 


1 .20E 


00 


6.30E- 


•01 


-7.60E- 


-01 


-1. 18E 


00 


-1.51E 


00 


-2.30E 


00 










1.81E 


00 


1.30E 


00 


7.90E- 


■01 


-6.40E-01 


-1.C7E 


00 


-1.42E 


00 


-2.24E 


00 










1.S7E 


.00 


1.44E 


00 


9.50E- 


■01 


-4.30E- 


-01 


-8.60E- 


-01 


-1.29E 


00 


-2.20E 


00 










-2.80E 


00 


-2.80E 


00 


-2.85E 


00 


-3.30E 


00 


-3.47E 


00 


-3.62E 


00 


-4. OOE 


00 










-2.40E 


00 


--2.46E 


00 


-2.65E 


00 


-3.10E 


00 


-3.27E 


00 


-3.41E 


00 


-3.81E 


00 










-1.84E 


00 


-1.88E 


00 


-2. OOE 


00 


-2.50E 


00 


-2.67E 


00 


-2.83E 


00 


-3.35E 


00 










-1.28E 


00 


-1.48E 


00 


-1.70E 


00 


-2.30E 


00 


-2.50E 


00 


-2. TOE 


00 


-3.24E 


00 










-1.90E- 


-01 


-3.30E- 


•01 


-5.30E- 


■01 


-1.41E 


00 


-1.71E 


00 


-2. OOE 


00 


-2.71E 


00 










1.17E 


00 


8.80E- 


•01 


4.30E-01 


-9.40E- 


-01 


-1.30E 


00 


-1.66E 


00 


-2.50E 


00 










1 .80E 


00 


1.43E 


00 


8.80E- 


■01 


-6.50E-01 


-1.C5E 


00 


-1.41E 


00 


-2.31E 


00 










2. lOE 


00 


1.61E 


00 


' I.OIE 


00 


-5.00E- 


-01 


-9.40E- 


-01 


-1.30E 


00 


-2.22E 


00 










2.69E 


00 


2.21E 


00 


1.55E 


00 


-2.50E- 


-01 


-7. OOE- 


•01 


-1.09E 


00 


-2. OOE 


00 











E 

00 

01 

01 

E 

01 

01 

E 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

02 

00 

00 

00 

01 

00 

GO 

01 

00 

00 

01 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

02 

00 

00 

00 

01 

00 

00 

01 

00 

00 

02 

00 

00 

01 

00 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 











TABLE 6 


CONT 


INUED 














1 


-2. 


. 51E 


00 


-2.36E 


00 


-2 


.25E 


00 


-2 


. 20E 


00 


-2 


.20E 


00 




-2. 


,25E 


00 


-2.30E 


00 


-2 


.37E 


00 


-2 


. A4E 


00 


-2 


.52E 


00 




-2. 


► 60E 


00 


-2.66E 


00 


-2 


.72E 


00 














1 


-1. 


, 84E 


00 


-1.62E 


00 


-1 


.54E 


00 


-1 


. 58E 


00 


-1 


.72E 


00 




-1. 


,89E 


00 


-2.03E 


00 


-2 


. 15E 


00 


-2 


.26E 


00 


-2 


.34E 


00 




-2. 


.41E 


00 


-2.48E 


00 


-2 


. 54E 


00 














1 


-1, 


,23E 


00 


-1.24E 


00 


-1 


.28E 


00 


-1 


. 36E 


00 


-i 


.49E 


00 




-1, 


,64E 


00 


-1.81E 


00 


-1 


.96E 


00 


-2 


.09E 


00 


-2 


.21E 


00 




-2. 


.31E 


00 


-2.40E 


00 


-2 


.47E 


00 














1 


-7, 


,90E- 


-01 


-8.20E- 


-01 


-8 


.90E- 


■01 


-1 


.C2E 


00 


-1 


.19E 


00 




-1. 


,38E 


00 


-1.53E 


00 


-1 


.69E 


00 


-1 


.82E 


00 


'1 


.96E 


00 




-2. 


,07E 


00 


-2.15E 


00 


-2 


.20E 


00 














1 


0. 


.0 




-I.IOE- 


-01 


-2 


.40E- 


•01 


-4 


. lOE- 


-01 


-6 


.lOE- 


-01 




-8. 


,30E- 


-01 


-1. 08E 


00 


-1 


.30E 


00 


-1 


.50E 


00 


-1 


.67E 


00 




-1. 


,79E 


00 


-1.89E 


00 


-1 


. 99E 


00 














1 


8. 


,20E- 


-01 


7.70E- 


-01 


5 


.80E- 


•01 


3 


.60E- 


-01 


2 


.OOE- 


-02 




-2. 


.90E- 


-01 


-5.90E- 


-01 


-8 


.50E- 


•01 


-1 


. lOE 


00 


-1 


.29E 


00 




-1. 


,44E 


00 


-1.54E 


00 


-1 


.63E 


00 














1 


1. 


.24E 


00 


l.llE 


00 


9 


.OGE- 


■01 


6 


. 50E- 


-01 


3 


.60E- 


-01 




1. 


,00E- 


-02 


-2.80E-01 


-5 


.60E- 


•01 


-8 


.50E- 


-01 


-1 


.llE 


00 




-1. 


,32E 


00 


-1.47E 


00 


-1 


.56E 


00 














1 


1. 


.87E 


00 


1.65E 


00 


1 


.35E 


00 


1 


.OlE 


00 


6 


.30E- 


-01 




2. 


.50E-01 


-I.IOE- 


-01 


-4 


.OOE- 


•01 


-6 


.80E- 


-01 


-8 


.90E- 


-01 




-1. 


,05E 


00 


-1.19E 


00 


-1 


.28E 


00 














1 


2. 


.50E 


00 


2.30E 


00 


1 


.85E 


00 


1 


.40E 


00 


8 


.90E- 


-01 




4. 


,40E- 


-01 


5.00E- 


-02 


-2 


.90E- 


•01 


-5 


. 80E- 


-01 


-7 


.70E- 


-01 




-9. 


.lOE-Ol 


-1.03E 


00 


-1 


.12E 


00 














1 


-1. 


.40E 


00 


-1.43E 


00 


-1 


.50E 


00 


-1 


.55E 


00 


-1 


.65E 


00 




-1. 


,70E 


00 


-1.78E 


00 


-1 


.82E 


00 


-1 


. 90E 


00 


-1 


.92E 


00 




-1. 


.98E 


00 


-2. OlE 


00 


-2 


.05E 


00 














1 


-1< 


.OlE 


00 


-1.15E 


00 


-1 


.25E 


00 


-1 


.35E 


00 


-1 


.45E 


00 




-1. 


.55E 


00 


-1.65E 


00 


-1 


.75E 


00 


-1 


.80E 


00 


-1 


.85E 


00 




-1, 


.86E 


00 


-1.90E 


00 


-1 


.95E 


00 














1 


-8. 


.50E- 


-01 


-9.00E- 


-01 


-9 


.50E-01 


-1 


.OOE 


00 


-1 


.lOE 


00 




-1. 


,25E 


00 


-1.35E 


00 


-1 


.45E 


00 


-1 


.55E 


00 


-1 


.60E 


00 




-1, 


.65E 


00 


-1.75E 


00 


-1 


.80E 


00 














1 


-5. 


.50E- 


■Oi 


-6.00E- 


-01 


-7 


.OOE- 


•01 


-8 


. 50E- 


-01 


-9 


.50E- 


-01 




-1, 


,10E 


00 


-1.20E 


00 


-1 


.30E 


00 


-1 


. 40E 


00 


-1 


.45E 


00 




-1, 


,55E 


00 


-1.65E 


00 


-1 


.78E 


00 














1 


0, 


.0 




-1.20E- 


-01 


-2 


.70E- 


•01 


-4 


.30E-01 


-5 


.80E- 


-01 




-7. 


,80E- 


-01 


-9.50E- 


-01 


-1 


.lOE 


00 


-1 


.22E 


00 


-1 


.36E 


00 




-1, 


, 46 E 


00 


-1.54E 


00 


-1 


.61E 


00 














1 


5. 


.lOE- 


-01 


4.20E- 


-01 


2 


.80E- 


■01 


8 


.OOE- 


-02 


-1 


.80E- 


-01 




-4, 


.50E- 


-01 


-6.20E- 
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-1 


.15E 


00 














1 


1 


-3. 


,00E- 


-02 


-5.00E- 


-02 


-9 


.OOE-02 


-1 


.70E- 


-01 


-2 


.20E-01 






-2, 


.50E-01 


-2.80E-01 


-3 


.80E- 


-01 


-5 


.OOE-01 


-6 


.30E- 


-01 






-7. 


.70E- 


-01 


-9.20E- 


-01 


-1 


.09E 


00 














1 


1 


0. 


.0 




-l.OOE- 


-02 


-5 


.OOE- 


-02 


-1 


.30E- 


-01 


-2 


.lOE-01 






^2, 


.80E- 


-01 


-2.90E- 


-01 


-5 


.OOE-01 


-6 


.20E- 


-01 


-7 


.60E- 


-01 






-9, 


,20E- 


-01 


-1.08E 


00 


-1 


.27E 


00 














1 


1 


6. 


,00E- 


-02 


4.00E-02 


1 


.OOE- 


-02 


-2 


.OOE- 


-02 


-1 


.lOE- 


-01 






-1, 


.80E- 


-01 


-2.80E- 


-01 


-2 


.90E- 


-01 


-5 


.40E- 


-01 


-6 


.50E- 


-01 






-7. 


. 50E- 


-01 


-8.80E-01 


-1 


.08E 


00 












-02 


1 


1 


8. 


,00E- 


-02 


5.00E- 


-02 


4 


.OOE- 


-02 


^1 


. OOE- 


-02 


-5 


.OOE- 






-1. 


,60E- 


-01 


-2.40E- 


-01 


-3 


.50E- 


-01 


-4 


. 50E- 


-01 


-5 


.80E- 


-01 






-7, 


,OOE- 


-01 


-8.20E- 


-01 


-9 


.50E- 


-01 














1 


1 


1. 


.30E- 


-01 


1.20E- 


-01 


8 


.OOE-02 


2 


.COE- 


-02 


-3 


.OOE- 


-02 






-1, 


.20E- 


-01 


-2.20E- 


-01 


-2 


.20E- 


-01 


-4 


. 30E- 


-01 


-5 


.20E- 


-01 






->6. 


,00E- 


-01 


-6.70E- 


-01 


-7 


.lOE- 


-01 














1 


1 


1. 


.30E- 


-01 


1.20E- 


-01 


1 


.OOE- 


-01 


5 


.COE- 


-02 


0 


.0 








-9. 


. OOE- 


-02 


-1.80E- 


-01 


-2 


.90E- 


-01 


-4 


. lOE- 


-01 


-5 


.OOE- 


-01 






-5. 


,70E- 


-01 


-6.00E- 


-01 


-5 


.80E- 


-01 
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TABLE 7: 



Data Set 1; 5445 points; Evaluation time = 15 



Program 


Type of rms 

Approximation error 


Maximum 

error 


Points 


Evaluation 

time 


LSTSQ 


Quad, 


,nd 

3 rd 


var. .005 


.029 


3267 


21 


LSTSQ 


Quad, 


var. .004 


.026 


1815 


21 


LSTSQ 


Quad, 


var. .004 


.016 


1485 


21 


LSTSQ 


Quad, 


4^^ 


var. .023 


.107 


1485 


21 


MINMAX 


Sum 1 


var. 


ftns 


1.379 


36 


4 


MINMAX 


Sum 2 


var. 


ftns 


.402 


474 


18 


APPROX 


Sum 1 


var. 


ftns .429 


2.582 


36 


4 


APPROX 


Sum 2 


var. 


ftns .113 


.669 


474 


18 


APPROX 


Sum 3 


var. 


ftns .058 


.148 


2684 


28 


lBLE 8 : 


Data set 


2; 1512 points; 


evaluation time = 


7 





Program 


Type of 
Approximation 


rms 

error 


Maximum 

error 


Points 


Evaluation 

time 


LSTSQ 


Linear, 1^^ var. 
Linear, 2^^ var. 
Linear, 3^^ var. 
Quad , 1 ^^ var . 


.110 


.416 


144 


6 


LSTSQ 


.065 


.249 


249 


6 


LSTSQ 


.079 


.230 


378 


6 


LSTSQ 


.023 


.126 


216 


9 


LSTSQ 


Quad, 2^^ var. 

Quad , 3 var . 


.028 


.128 


504 


9 


LSTSQ 


.052 


.262 


567 


9 


MINMAX 


Sum 1 var. ftns. 




.370 


38 


3 


MINMAX 


Sum 2 var. ftns. 




.156 


429 


9 


APPROX. 


Sum 1 var. ftns. 


.137 


.509 


38 


3 


APPROX. 


Sum 2 var. ftns. 


.073 


.252 


429 


9 



TABLE 9; Data Set 3; 1512 points; evaluation time ■ 7 



Program 


Type of 
Approximation 


rms 

error 


Maximum 

error 


Points 


Evaluation 

time 


LSTSQ 


Quad, 1^^ var. 
Quad, 2^^ var. 


.028 


.160 


216 


9 


LSTSQ 


.027 


.094 


504 


9 


LSTSQ 


Quad, 3 var. 


.035 


.140 


567 


9 


MINMAX 


Sum 1 Var. ftns 




.448 


38 


3 


MINMAX 


Stun 1 Var. ftns 




.164 


429 


9 


APPROX 


Sum 1 Var. ftns 


.137 


.523 


38 


3 


APPROX 


Sum 2 Var. ftns 


.075 


.219 


429 


9 



TABLE 10: 


Data Set 4; 1512 ] 


points; 


evaluation time = 


7 




Program 


Type of 


rma 


Maximum 


Points 


Evaluation 




Approximation 


error 


error 




time 


LSTSQ 


Quad, 1^^ Var. 


.093 


.519 


216 


9 


LSTSQ 


Quad, 2^, Var. 


.049 


.285 


504 


9 


LSTSQ 


Quad, 3^ Var. 


.101 


.375 


567 


9 


MINMAX 


Sum 1 Var. ftns 




1.363 


38 


3 


MINMAX 


Sum 2 Var. ftns 




.514 


429 


9 


APPROX 


Sum 1 Var. ftns 


.318 


2.581 


38 


3 


APPROX 


Sum 2 Var. ftns 


.133 


1.083 


429 


9 



I4QJ 



TABLE 11: 
Program 



LSTSQ 

LSTSQ 

LSTSQ 

MINMAX 

MINMAX 

APPROX 

APPROX 



TABLE 12: 



LSTSQ 

LSTSQ 

LSTSQ 

MINMAX 

MINMAX 

APPROX 

APPROX 



Data Set 5; 1512 points; evaluation time = 7 



Type of 


rms 


Maximum 


Points 


Evaluation 


Approximation 


error 


error 




time 


Quad, 


2d Var. 
3'^'* Var. 


.079 


.671 


216 


9 


Quad, 


.054 


.237 


504 


9 


Quad, 


.112 


.346 


567 


9 


Sum 1 


Var. ftns. 




1.515 


38 


3 


Sum 2 


Var. ftns. 




.988 


429 


9 


Sum 1 


Var. ftns. 


.442 


2.687 


38 


3 


Sum 2 


Var. ftns. 


.341 


1.248 


429 


9 


Data set 


6; 1170 points; 


evaluation time = 


7 





Quad, 


^st 

1 , Var. 
2”'! Var. 
3''^ var. 


.069 


.448 


270 


Quad, 


.066 


.398 


390 


Quad, 


.166 


.850 


351 


Sum 1 


Var. ftns 




1.623 


32 


Sum 2 


Var. ftns 




.549 


337 


Sum 1 


Var. ftns 


.567 


2.179 


32 


Sum 2 


Var. ftns 


.210 


.802 


337 
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A. 0 Program Descriptions 



In this appendix the usage of programs MINMAX and LSTSQ are described, 
including input and output description, and a sample run for each. 

A. 1 Description of Program MINMAX 

Program MINMAX is a Fortran IV version of the algorithm outlined in 
Section 3.1. The program package consists of the main program and two sub- 
routines. Subroutine RDN(N) is supplied by the user and is used to read or 
generate the input function (table). The input argument N is the array of 
the number of entries in each variable. The version of RDN given in the 
program listings is the one used to read the test data supplied by NMC. 
Subroutine FAP is used to find the maximum and minimum values of an array, 
over specified limits of the subscripts. The average of these values are 
subtracted from each member of the array searched. 

The common statement 

COMMON F(l,2,3,4,5) 

must be supplied by the user and appears in the main program and both sub- 
routines. The actual dimensions used must be at least as large as the array 
to be input. 

The input to the main program consists of one card punched in format 
(715), The description of the input values follows. A blank card ends 



the computer run. 
Columns 

(right adjusted) 


Program 

Symbol 


Variable 




Comments 


1-5 


N(l) 




Number of 


entries for 1st Variable 


6-10 


N(2) 




Number of 


entries for 2nd Variable 


11-15 


N(3) 




Number of 


entries for 3rd Variable 


16-20 


N(4) 




Number of 


entries for 4th Variable 


21-25 


N(5) 


A 


Number of 


entries for 5th Variable 


26-30 


NPR 




NPR ^ 0 will cause printout of the 
details of each iteration. 

Normally NPR = 0 


31-35 


NVF 




Number of 
functions ; 


variables for the fitting 
1 or 2 



MAIN PROGRAM 



M I N M A X 



A. 2 

MAINPROGRAM MINMAX 

DIMfcNSICN F1(5,25J ,N(5) ,M( 5J ,MX(5) »F2( 10, 25,251 

EQUIVALENCE ( N( 1 ) ,N 1 ) , ( N ( 2 ) , N2 ) , ( N ( 3 ) ,N3 ) , ( N (^ ) , N A ) , ( N ( 5 ) , N5 ) 

COMMON F(26,10, 10,1, 1) 



100 



150 



160 

170 



190 

200 



210 

225 



PROGRAM MINMAX OBTAINS A SUM OF 
FUNCTION OF UP TO FIVE VARIABLES 



FUNCTIONS APPROXIMATION TO A 
WHICH MINIMIZES THE MAXIMUM 



SUBROUTINE RDNIN) IS USER SUPPLIED 
UR READ THE DATA IN THE ARRAY F 



AND IS USED TO GENERATE 



THE USER MUST ALSO SUPPLY APPROPRIATE STATEMENT OF THE FORM 
COMMON F(l,2,3,4,5) 

WHERE THE ACTUAL DIMENSIONS USED ARE ADEQUATE FOR THE INPUT 
THIS COMMON STATEMENT MUST APPEAR IN PROGRAM MINMAX AS WELL 
SUBROUTINES RDN AND FAP 



F 

AS 



IN 



INPUT VARIABLE 
(N( I ) , I =1,5J = 
FIRST, SECOND, 



DESCRIPTION, ONE (1) CARD, FORMAT! 715) 

N1 ,N2,N3,N4,N5 ARE THE NUMBER OF ENTRIES 
ETC VARIABLES FOR THE INPUT FUNCTION F. 



IN The 



NPR = 
NPR = 
ERROR 

NVF = 
NVF = 



0 IS NORMAL 

1 WILL GIVE A PRINTOUT 
AT EACH ITERATION 



OF THE FITTING FUNCTIONS AND THE 



1 GIVES AN APPROXIMATION BY A SUM OF ONE VARIABLE FUNCTIONS 

2 GIVES AN APPROXIMATION BY A SUM OF TWO VARIABLE FUNCTIONS 



250 

300 



READ(5,1)N,NPR,NVF 
IF(Nl.LE.O) STOP 
CALL RDN(N) 

WRITE(6,5) 

DO 150 15=1, N5 
DO 150 14=1, N4 
DO 150 13=1, N3 
DC 150 12=1, N2 

WRITE(6,3)I2,I3,I4,I5,(F(I1,I2,I3,I4,I5),I1=1,N1) 

00 160 1=1,5 

IF(N( I ) .EQ. DGO TO 170 

NV = I 

NVM = NV - 1 
ITMAX = 20 
IT = 1 
EPS = .5E-4 
DO 200 1=1,5 
DC 190 J=l,25 
FI (I , J) = 0.0 
M( I ) = 1 
MX(I) = NdJ 
DC 210 1=1,10 
DO 210 J=l,25 
DO 210 K=l,25 
F2(I,J,K) = 0. 

DX = 0.0 

IF(NVF.GT.l )GU TO 310 

THE CALCULATION OF ONE VARIABLE APPROXIMATIONS 

DO 300 I=1,NV 
NP = N( I ) 

DC 250 J=1,NP 
M( I ) = J 
Mxn ) = J 
CALL FAP(FV,M,MX) 

DX = AMAXK ABS(FV) ,DX) 

F1(I,J) = FV + F1(I,J) 

M( I) = 1 
MX(I) = NP 
CONTINUE 
GO TO 390 



MINOOOlO 
MIN00020 
MIN00030 
MIN00040 
MIN00050 
MIN00060 
ERR0RMIN00070 
MIN00080 
MIN00090 
MINOOlOO 
MINOOllO 
MIN00120 
MIN00130 
MIN00140 
MIN00150 
MIN00160 
MINC0170 
MIN00180 
KIN00190 
MIN00200 
MIN00210 
MIN00220 
MIN00230 
MIN00240 
MIN00250 
MIN00260 
KIN00270 
MIN00280 
MIN00290 
MIN00300 
MIN00310 
MIN00320 
MIN00230 
MIN00340 
MIN00350 
MIN00360 
MIN00370 
MIN00380 
MIN00390 
MIN00400 
MIN00410 
MIN00420 
MIN00430 
MIN00440 
MIN00^50 
MINC0460 
MIN00470 
MIN00480 
MIN00490 
MIN00500 
MINC0510 
MIN00520 
MIN00530 
MIN00540 
MIN00550 
MIN00560 
MIN00570 
MIN00580 
MIN00590 
MIN00600 
MIN00610 
MINC0620 
MIN00630 
MIN00640 
MIN00650 
MIN00660 
MINC0670 
MIN00680 
MIN00690 
MIN00700 
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THE CALCULATION OF TWO VARIABLE APPROXIMATIONS 



310 MF = 0 

DC 370 1=1, NVM 
IX = I + 1 
NPI = N(I) 

DO 360 K=IX,NV 
MF = MF + 1 
NPK = N(K) 

Du 350 JI= 1,NPI 
M(I) = JI 
MXd ) = JI 
DO 340 JK=1,NPK 
M(K» = JK 
MX(K) = JK 
CALL FAP(FV,M.MXJ 
DX = AMAXK ABS(FVJ,DX) 

340 F2(MF,JI,JK) = F2(MF,JI,JK) + FV 
350 CONTINUE 
M(K) = 1 
MX(K) = NPK 
360 CONTINUE 
M( I ) = 1 
MXd) = NPI 
370 CONTINUE 
390 CONTINUE 

IF(DX.LT.EPS)GO TO 400 
IT = IT + 1 

IF (IT.GT.ITMAX)GO TO 500 
IF(NPR.EC.O)GO TO 225 
400 WRITE(6,7)DX, IT 

IFINVF.GT.l )G0 TO 412 
WRITE(6,6) 

DO 410 1=1, NV 
NP = N( I J 

410 WRIT£(6,8) I ,(F1( I , J) ,J = 1,NP) 

GO TO 418 
412 MF = 0 

WRITE(6,6) 

DC 415 1=1, NVM 
NPI = Nd ) 

IX = I + 1 
DC 415 K=IX,NV 
WRITE(6,12)K 
MF = MF + 1 
NPK = N(K) 

DO 415 JK=1,NPK 

415 WRITE(6,9)I,K,JK, (F2(MF,IK, JK ),IK=1,N0I) 

418 CONTINUE 
AM = 0. 

DO 419 15=1, N5 
DO 419 14=1, N4 
DO 419 13=1, N3 
DO 419 12 = 1, N2 
DO 419 11=1, N1 

419 AM = AMAXK ABS(F( II, 12, 13, 14, 15) ) ,AM) 

WPITE(6.11) AM 
DO 420 15=1, N5 
DO 420 14=1, N4 
DO 420 13=1, N3 

■ DO 420 12=1, N2 

420 WRITE(o,3)I2,I3,l4,I5,(Fdl,I2,I3,I4,I5),Il = l,Nl) 
IF(0X.LT.EPS.0R.IT.GT.ITMAX)G0 to 100 
GO TO 225 

500 WRITE(6,4)IT 
GO TO 400 

1 FGPMAT(16I5) 

3 FCRMAT(4I4,lPlQE11.3/i 16X , lOE 1 1 . 3 ) ) 

4 F0RMAT(45H0C0NVERGENC£ FAILURE, NUMBER OF ITERATIONS IS ,15) 

5 F0RMAT(64H1THE INPUT FUNCTION, FIRST INDEPENDENT VARIATION ACROSS 

ITHE PAGE/16H0 12 13 14 15) 

6 FORMAT (91-OFUNCTION) 

7 FORMAT! 75H0THE APPROXIMATING FUNCTIONS FOLLOW. LAST CF AND NUM3E RM I NO 1440 

1 OF ITERATIONS WERE , IPE 12 . 5 , 1 10// ) MIN01450 



MIN00710 

MIN00720 

MIN00730 

MIN00740 

MIN00750 

MIN00760 

MIN00770 

MINQ0780 

MIN00790 

MIN00800 

MIN00810 

MIN00820 

MIN00830 

KIN00840 

MIN00850 

MIN00360 

MIN00870 

MIN00880 

MIN00890 

MIN00900 

MINC0910 

MIN00920 

MINC0930 

MIN00940 

MIN00950 

MINC0960 

MIN00970 

MIN00980 

MINC0990 

MINOIOOO 

MINOlOlO 

MIN01020 

MIN01030 

MIN01040 

MIN01050 

MIN01060 

MIN01070 

MIN01080 

MIN01090 

MINOllOO 

MINOlllO 

MIN01120 

MIN01130 

MIN01140 

MIN01150 

MIN01160 

MIN01170 

MIN01180 

MIN01190 

MIN01200 

MIN01210 

MIN01220 

MIN01230 

MIN01240 

MIN01250 

MIN01260 

MIN01270 

MIN01280 

MIN01290 

MIN01300 

MIN01310 

MIN01320 

MIN01330 

MIN01340 

MIN01350 

MIN01360 

MIN01370 

MIN01380 

MIN01390 

MIN01400 

MIN01410 

MIN01420 

MING1430 
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8 FCRMAT( 2HOF,Il,7X,lP10E11.3/(10X,10E11.3n MIN01460 

9 F0RMAT(2h F ,2 ll . 16, IPlOEll . 3/UOX, lOE 11. 3) ) MINC1470 

11 FCRMAK A3Hlfh£ ERROR ARRAY FCLLQ^JS. MAXIMUM ERROR IS,1PE12.4/ MIN01480 

1 16H0 12 13 14 15) MIN01490 

12 FORMAT (/8X, IFI, II ) MIN01500 

END MIN01510 



SUBROUTINE RDN 

SUBROUTINE RON(N) 

DIMENSION N(l) 

COMMON F(26,10,10,l, 1) 

N1 = N( 1 ) 

N2 = N(2) 

N3 = N ( 3 ) 

N4 = N(4) 

N5 = N(5) 

DO 110 15=1, N5 
DO 110 I<i=l,N4 
DO 110 13=1, N3 
DC 110 12=1, N2 

110 READ(5,2)(F(I1,I2,I3,!4,I5),I1=1,N1) 

RETURN 

2 FCRMAT(9( 1X,F5.2,2XJ ) 

END 



RDNOOOlO 

RDN0002C 

RDN00030 

RDN00040 

RDN00050 

RDN00060 

RDN00070 

RDN00080 

RDN00090 

RDNOOlOO 

RDNCOllO 

RDN00120 

RDN00130 

RDMG0140 

RDN00150 

RDN00160 



SUBROUTINE FAR 

SUBROUTINE FAP(FV,M,MX) 

DIMENSION M(l),MXm 
COMMON F(26, 10, 10,1,1) 

FAP FINDS THE MAX AND MIN VALUES OF F WITH INDEX II VARYING 
FROM M(l) TO MX(1), 12 VARYING BETWEEN M(2) AND MX(2), ETC. 

FV IS THE AVERAGE VALUE OF THE MAXIMUM AND MINIMUM VALUES , AND 
SUBTRACTED FROM THAT SET OF FUNCTION VALUES 



IBl = M(l) 

162 = M(2) 

IB3 = M(3) 

IB4 = M(4) 

IBS = M<5) 

N1 = MX(1) 

N2 = MX(2) 

N3 = MX 13) 

N4 = MX(4) 

N5 = MX(5) 

A = -I.EIO 

Z = l.ElO 

CO 100 I1=IB1,N1 

CO 100 12=IB2,N2 

DC 100 I2=IB3,N3 

DC 100 I4=IB4,N4 

DC 100 15 = IB5,N5 

A = AMAXHFdl, 12, 13,14,15) , A) 

100 Z = AMIM(F(I1,I2,I3,I4,I5),Z) 

FV = (A + Z)/2. 

DO 120 I1=IB1,N1 
DC 120 I2=IB2,N2 
CO 120 I3=IE3,N3 
DC 120 IA=I,B4,N4 
DO 120 I5=IB5,N5 

120 F( 11,12, 13, 14,15) = F(Il,I2,I3,d,I5) 
RETURN 
END 



FAPOOOlO 
FAP00020 
FAP00030 
FAP00040 
FAP00050 
FAP00060 
FAP00C70 
FAP00080 
FAP00090 
FAPOO 100 
FAPOOllO 
FAP00120 
FAP00130 
FAP00140 
FAP00150 
FAP00160 
FAPOOITO 
FAP00180 
FAP00190 
FAP00200 
FAP00210 
FAP00220 
FAP00230 
FAP00240 
FAP00250 
FAP00260 
FAP00270 
FAP00280 
FAP00290 
FAPC0300 
FAP00310 
FAP00320 
FAP00330 
FAP00340 

- FV FAP00350 

FAP00360 

FAP00370 
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A. 3 Sample Data and Output for MINMAX 



The following data set was run on the program as listed, a*.d the 
output follows. It should be noted that this data is for illustrative 
purpose only. 



(Input data for MINMAX, as punched ) 



4 3 2 

1.0 2.0 

-.5 1.5 

^ 0 1.0 

. 0.0 

. . 5 

.5 ^.2 

4 3 2 

0 2.0 

5 1.5 

0 1.0 

. 0.0 

. -.5 

.5 -.2 



110 

3.5 4.0 

2.9 3.8 

2.2 2.5 

1.5 2.5 

0.5 2.0 

0.3 1.9 

110 

3.5 4.0 

2.9 3.8 

2.2 2.5 

1.5 2.5 

0.5 2.0 

0.3 1.9 



1 



2 



(Output for MINMAX) 

THE INPUT FUNCTION, FIRST INDEPENDENT VARIATION ACROSS THE PAGE 



12 


13 


14 


15 














1 


1 


1 


1 


l.OOOE 


00 


0.0 


5.000E-01 


0.0 




2 


1 


1 


1 


-5.000E- 


•01 


1.500E 00 


2.900E 00 


3.800E 


00 


3 


1 


1 


1 


0.0 




l.OOOE 00 


2.200E 00 


2.500E 


00 


1 


2 


1 


1 


-l.OOOE 


00 


0.0 


1.500E 00 


2.500E 


00 


2 


2 


1 


1 


-2.000E 


00 


-5.000E-01 


5.000E-01 


2.000E 


00 


3 


2 


1 


1 


-1.500E 


00 


2.000E-01 


3.000E-01 


9.000E- 


-01 



THE APPROXIMATING FUNCTIONS FOLLOW. LAST DF AND NUMBER OF-^ 



FUNCTION 

FI -2.000E-01 4.500E-01 1.650E 00 1.900E 00 

F2 -3.500E-01 3.500E-01 -7.501E-02 

F3 0.0 -6.000E-01 




ITERATIONS WtRE 2.06232E-05 
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THE ERROR ARRAY FOLLOWS. MAXIMUM ERROR IS 1.5500E 00 



12 13 I^ 

1 1 1 

2 11 
3 11 

12 1 
2 2 1 

3 2 1 



15 

1 

1 

1 

1 

1 

1 



1.550E CO 
-6.500 E-01 
2.750E-01 
1.500E-01 
-1.550E 00 
-6. 250E-01 



-l.OOOE-01 
7.000E-01 
6.250E-01 
5.000E-01 
-7.000E-01 
4.25 0E-01 



-8.000E-01 

9.000E-01 

6.250E-01 

8.000E-01 

-9.000E-01 

-6.750E-01 



-1.550E 00 
1.550E 00 
6. 750 E-01 
1.550E 00 
3.500E-01 
-3.250E-01 



THE INPUT FUNCTION, FIRST INDEPENDENT VARIATION ACROSS THE PAGE 



12 

1 

2 

3 

1 

2 

3 



13 

1 

1 

1 

2 

2 

2 



lA 

1 

1 

1 

1 

1 

1 



15 

1 



0.0 5.000E-01 0.0 

1.500E 00 2.900E 00 3.800E 00 

l.OOOE 00 2.200E 00 2.500E 00 

0.0 1.500E 00 2.500E 00 

-2.000E 00 -5.000E-01 5.000E-01 2.000E 00 

-1.500E 00 2.000E-01 3.000E-01 9.000E-01 



l.OOOE 00 
1 -5.000E-01 
1 0.0 

1 -l.OOOE 00 
1 
1 



THE APPROXIMATING FUNCTIONS FOLLOW. LAST DF AND NUMBER OF- 



FUNCTION 






ITERATIONS WERE 0.0 



12 



F12 


1 


0.0 


0.0 


l.OOOE 00 


1.250E 


00 


F12 


2 


-1.250E 00 


5.000E-01 


1.700E 00 


2.900E 


00 


F12 


3 


-7.500E-01 


6.000E-01 


1.250E 00 


1.700E 


00 


F13 


13 

1 


8.750E-01 


2.500E-01 


3.500E-01 


-1.750E 


-01 


F13 


2 


-8.750E-01 


-2.500E-01 


-3 .500E-01 


1.750E 


-01 


F23 


13 

1 


-4.750E-01 


4.750E-01 


4.250E-01 






F23 


2 


4.750E-01 


-4.750E-01 


-4.250E-01 







THE 


ERROR 


ARRAY FOLLOWS. MAXIMUM ERROR IS 6.0000E-01 


12 


13 


14 


15 






1 


1 


1 


1 6.000E-01 2.250E-01 


-3.750E-01 


-6.000E-01 


2 


1 


1 


1 -6.000E-01 2.750E-01 


3.750E-01 


6. OOOE-01 


3 


1 


1 


1 -5.500E-01 -2.750E-01 


1.750E-01 


5. 500E-01 


1 


2 


1 


1 -6.000E-01 -2.250E-01 


3.750E-01 


6.000E-01 


2 


2 


1 


1 6.000E-01 -2.750E-01 


-3.750E-01 


-6. OOOE'Ol 


3 


2 


1 


1 5.500E-01 2.750E-01 


-1.750E-01 


-5.500E-01 



The output is mostly self explanatory. Output of the input function and 
error arrays are in the same format as the information in Tables 6. 

The *'LAST DF" is the value of max |Af| during the last iteration (last 
value of ty as outlined in Section 3.1). The stopping criterion is that 
this value be less than EPS = 5*10 ^ . The maximum number of iterations 
allowed is ITMAX = 20. Both of these values are easily changed. 

A. 4 Description of Program LSTSQ 

Program LSTSQ is a Fortran IV version of the idea in Section 3.4. 

The program package consists of the main program, subroutines LLSQFT, 

WLSQ, RD, and function PHI. The latter two are supplied by the user. 

RD is used to read the variable and function data. The user specifies 
in Function PHI(I,X) up to three basis functions (J)^(x), and 

({J)^,or (j)^ and (j)^ may be identically zero to obtain approximations with 2 
or 1 basis function, respectively.) The program will then perform a weighted 
(weights user specified) least squares fit in one variable. This variable 
may be specified; a search may be made for the variable in which the fit is 
best, and the approximation computed in that variable; or the rms and maximum 

errors may be computed for each variable. 

Subroutine LLSQFT picks the appropriate data from the arrays F,X, and W, 
and then calls subroutine WLSQ which does a weighted least squares fit on 
the data. This is done by first computing an orthogonal basis, and then 
transforming back to the basis ^2.'^2*^3 ’ 

The user must also supply the proper common statement of the form 
COMMON/TABLE/Nl,N2,N3,N4,N5,N6,F(l,2,3,4,5,6) ,X(25,6) ,W(25,6) 

The actual dimensions of F must be sufficient to store the input F array, 
plus 3 more than that for any variable in which the approximation is to be 
output, since the coefficients of ())^, <i>^y and (j)^ are stored in these extra 
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locations. The second diiuension of X and W njust he as large as 

Max N. . The common statement appears in the main program and sub- 
l£i<_6 ^ 

routines RD and LLSQFT. 

The input to LSTSQ consists of a variable number of cards, depend- 
ing on the option concerning user specified weights. The first card is 
punched in format (1315), as follows. A blank card ends the computer run. 



Columns Program Variable 

(right adjusted) Symbol 



Comments 



1-5 


N(l) 


6-10 


N(2) 


11-15 


N(3) 


16-20 


N(4) 


21-25 


N(5) 


26-30 


N(6) 


31-35 


NW(1) 


36-40 


NW(2) 


41-45 


NW(3) 


46-50 


NW(4) 


51-55 


NW(5) 


56-60 


NW(6) 


61-65 


NC0DE 



N, 

Nc 

n: 



Number of entries for 1st variable 
Number of entries for 2nd variable 
Number of entries for 3rd variable 
Number of entries for 4th variable 
Number of entries for 5th variable 
Number of entries for 6th variable 
^NW(I) = 0 indicates weights are to 
I be all equal to one. NW(I) ^ 0 
'indicates weights are to be input 
for the varible > as card two. 



/n( 



NC0DE = K (1 < K £ 6) indicates 
\a fit on the variable is to be 
■performed. NC0DE = 0 indicates a 
y search for rms and maximum errors 
y for the fit in each variable. 

I NC0DE < 0 indicates a search as 
I above, and the approximation output 
y for smallest rms error. 



If the weights are to be input, they are read in order under format (8E10.0). 



More than one card may be necessary for each set of weights read. Each set 
begins on a new card. 
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A. 5 



MAINPROGRAM 



L S T S Q 



MAINPROGRAM LSTSO 

COMMON/ TABL E/Nl, N2, N3, N4, N5 ,N6, F( 25, 12,il,i,l,l»,X(25,6),W(25,6) 
DIMENSION N(6),NW(6) 

EQUIVALENCE ( N( 1 ) , Ni ) , ( Nw ( 1 ) , K1 ) , ( NN( 2 ) ,K 2 ) , (NW( 3 ) , K3 ) , I NW ( A ) 

1 (NW(5) ,K5) ,(NW(6) ,K6) 

PROGRAM LSTSG PERFORMS A w/EIGHTED LEAST SQUARES CURVE FIT IN L 






VARIABLE OF 
INDEPENDENT 



MULTI DIMENSIONAL 
VARIABLES. 



TABLE. THERE MAY BE UP TO SIX 



SUBROUTINE RD MUST BE SUPPLIED BY THE USER 
INDEPENDENT AND DEPENDENT VARIABLE ARRAYS. 

A STATEMENT OF THE FORM 



AND READS THE 



COMMON /TAdL£/Nl,N2,N3,N4,N5,N6,F(l,2,3,4,5,6) , X( M , 6 ) , is ( M , 6 ) 

MUST BE SUPPLIED BY THE USER. THE DIMENSIONS OF F MUST BE 
ADEQUATE FOR THE INPUT F. THE FIRST DIMENSION OF X AND W MUST 
BE AS LARGE AS THE LARGEST DIMENSION OF F. THE COMMON STATEMENT 
APPEARS IN SUBROUTINES RO AND LLSQFT AS WELL AS THE MAIN PROGRAM. 
X( J, I ) , J = i,NI AND W(J,n,J = l,NI ARE THE ABSCISSAS AND WEIGHTS 
FOR THE ITH VARIABLE. 

THE FITTING FUNCTIONS MUST BE SUPPLIED IN THE FORM OF A FUNCTION 
TYPE SUBROUTINE. IF ONLY ONE OR TWO FUNCTIONS ARE TO BE USED IN 
THE APPROXIMATION, THEN PHI(2,X) AND PHI(3,X) OR PHI(3,X), 
RESPECTIVELY, ARE SET TO ZERO FOR ALL X. 



INPUT DESCRIPTION 
CARO ONE FORMAT! 1315) 

N1 ,N2,N3,N4,M5,N6 THE DIMENSIONS 
NW(I),I=1,6 NW(I) = 1 INDICATES 
IS TO BE INPUT. IF NW(I) = 0, THE 
NCODE (COLUMN 61-65) = -1 



OF THE INPUT ARRAY F 
THAT AN ARRAY OF N(I) WEIGHTS 
WEIGHTS ARE ALL SET t“ 



ARE ALL SET TO 1 
INDICATES A SEARCH IS TO BE MADE 
FOR THE VARIABLE THAT YIELDS THE 
APPROXIMATION WITH THE SMALLEST 
ERROR, AND THE APPROXIMATION TO 
OUTPUT. 



LSTOOOlO 
LST00020 
LST00030 
, LST00040 
LST00050 
LST00060 
LST00070 
LST00080 
LST00090 
LSTOOlOO 
LSTOOllO 
LST00120 
LST00130 
LST00140 
LST00150 
LST00160 
LST00170 
LST00180 
LST00190 
LST00200 
LST00210 
LST00220 
LST00230 
LST00240 
LST00250 
LST00260 
LST00270 
LST00280 
LST00290 
LST00300 
LST00310 
LST00320 
LST00330 
LST00340 
LST00350 
RMSLST00360 
BE LST00370 
LST00380 



= 0 INDICATES THAT THE APPROXIMATION I SLST00390 

TO BE COMPUTED FDR EACH VARIABLE LST00400 
AND THE RMS AND MAXIMUM ERROR GIVENLST00410 
= K INDICATES THAT AN APPROXIMATION IN LST00420 

THE KTH VARIABLE IS TO BE COMPUTED 



CARD 

THIS 

CARD 

THAN 

THE 



TWO F0RMAT(8E10.0 ) 

CARD IS REPEATED ONE FOR EACH 
MAY APPEAR ZERO TO SIX TIMES. 
ONE CARD IF N( I) IS 
VALUES OF THE WEIGHT 



NON-ZERO NW(I). THUS. THIS 
EACH 'CARD* MAY CONSIST OF 



WHICH NW(I ) IS NON-ZERO 



100 READ(5,1)N,NW, NCODE 
IF(Nl.LE.O) STOP 
WRITE(6,7) 

M = MAXO(N( 1) ,N(2) ,N(3) ,N(4) ,N(5) ,N(6 )) 
DO 105 J=l,e 
DO 105 1=1, M 
105 W( I , J) = 1. 

DO 110 J=l,c 

IF(NW( J ) .EQ.O)GU TO 110 

M = N(J) 

READ(5,4)(W(I,J) ,I=1,M) 

110 CONTINUE 
CALL RD 
DC 120 J=l,o 
M = N(J) 

IF(M.EQ.1)GQ TO 120 

WRITE(6,9)J,(W(I,J),I=1,M) 

WRITE(6,5)J,{X(I,J),I=1,M) 



LST00430 
LST00440 
LST00450 
LST00460 
M0RELST00470 

GREATER THAN EIGHT. LST00480 

ARE LISTED IN ORDER FOR THE VARIABLES FORLST00490 

LST00500 
LST00510 
LST00520 
LST00530 
LST00540 
LST00550 
LST00560 
LST00570 
LST00580 
LST00590 
LST00600 
LSTOOblO 
LST00620 
LST00630 
LST00640 
LST00650 
LST0(i660 
LST00670 
LST00680 
LST00690 
LST00700 
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120 



CONTINUE 

WRITE(6,8] 

DO 200 16=1, N6 



200 



DO 

DO 

DO 

DO 



200 

200 

200 

200 



300 

400 



500 

510 



520 



15=1, N5 
14=1, N4 
13=1, N3 
. 12=1, N2 

WRITE(6, 3)12, 13, 14, 15, 16, (Fill, 12, 13, 14, 15, 16), 11 = 1, NU 

IF(NCCOE.GT.O)GC TO 400 

RMSL = I.EIO 

DO 300 1=1,6 

IF(N( I) .LE.DGO TO 300 

CALL LLSCFTI I ,RMS ,ERM,0) 

WRITE(6,6) I ,RMS,tRM 
IFIRMS.GT.RMSDGO TO 300 
RMSL = RMS 
LR = I 
CONTINUE 

IF(NCODE.EQ.O)GO TO 100 
NCOOE = LR 

CALL LLSCFT(NC0DE,RMS,ERM,1J 
WRITE(6,11) 

DO 500 16=1, N6 
15=1, N5 
14=1 ,N4 
13 = 1, N3 
. 12 = 1 ,N2 

kNRITE(6,3)I2,I3,I4,I5,I6,(F(Il,I2,I3,I4,I5,I6),Il = l,Nl) 
DO 510 1=1,6 
Nwn ) = 1 

NW(NCODE) = N(NCuOE) + 1 
N(NCODE) = NW(NCOOE) + 2 

WRITEI6, 10) NCCDE, RMS, ERM, NCOOE, NW(NCODE) ,N (NCOOE) ,NCODE 
DO 520 I6=K6,N6 
I5=K5,N5 
I4=K4,N4 
I3=K3,N3 
I2=K2,N2 

WRITE(6,3)I2,I3,I4,I5,I6,(F(I1,I2,I3,I4,I5,I6),I1 = K1,M) 
GO TO 100 



DC 

DO 

on 

DO 



500 

500 

500 

500 



DO 

DO 

DO 

DO 



520 

520 

520 

520 



9 

10 



11 



FORMAT! 1615) 

FORMAT! 8E10.0) 

FORMAT ! 513, IP 10E 11.3/! 15X, 10E11.3)) 

FORMAT! 8E10.0) 

FORMAT !3H X , I 1 , 5X, IPlOE 1 1 . 3/ ! 9X , lOE 1 1 . 3 ) ) 

F0RMAT!8H0F0R TH£,I3,43HTH VARIABLE THE RMS AND MAXIMUM ERRORS 
1 ,1P2E13.4) 

F0RMAT!26H1THE INDEPENDENT VARIABLES//) 

F0RMAT(68H0THE DEPENDENT VARIABLE, FIRST INDEPENDENT VARIATION 
loss THE PAGE /15H0I2 13 14 15 16) 

FORMAT !3h WT , I 1 , 5X, 1 PlOEl 1 . 3/ ! 9X , lOE 1 1 . 3 ) ) 

FORMAT! 15H1THE FIT ON THE,I3,29HTH VARIABLE, WITH RMS ERROR =, 
1 1PE11.3 /32X,15HMAXIMUM ERROR = , E 1 1 . 3/58H0THE COEFFICIENTS 
2HI1, PHI2, AND PHIS ARE FUNCTIONS 0F/12H ALL BUT THE,I3,39HTH 
3A6LE, AND ARE LISTED HERE AS THE , 13 , 2HTH/ 12H THROUGH THE, 13, 

4 15HTH VALUES FOR I.I1,15H, RESPECTI VELY. / 1 5H0I 2 13 14 15 16) 

F0RMAT!24hlTHE ERROR ARRAY FOLLOWS/ 15H0I2 13 14 15 16) 

END 



LST00710 
LST00720 
LST00730 
LST00740 
LST00750 
LST00760 
LST00770 
LST00780 
LST00790 
LST00800 
LST00810 
LST00820 
LST00830 
LST00840 
LST00850 
LST00860 
LST00870 
LST00880 
LST00890 
LST00900 
LST00910 
LST00920 
LSTC0930 
LST00940 
LST00950 
LST00960 
LST00970 
LST00980 
LST00990 
LSTOIOOO 
LSTOlOlO 
LST01020 
LST01030 
LST01040 
LST01050 
LST01060 
LST01070 
LST01080 
LST01090 
LSTOl 100 
LSTOl 110 
LST01120 
LST01130 
LST01140 
LST01150 
LSTC1160 
ARELST01170 
LST01180 
LSTOl 190 
ACRLST01200 
LST01210 
LST01220 
LST01230 
OF PLST01240 
VARILST01250 
LST01260 
LST01270 
LST01280 
LST01290 



SUBROUTINE RD 

SUBROUTINE RO 

CCMM0N/TABLE/N1,N2,N3,N4,N5,N6, F!25, 12,11,1,1,1),X!25,6),W!25,6) 
DO 100 13=1, N3 
DO 100 12=1, N2 

100 REA0!5,1)!F(I1,I2,I3,1,1,1),X!I1,1),I1=1,9),X!I3,3),X!I2,2), 

1 !F!I1, I2,I3,1,1,1),X!I1,1),I1=10,N1) 

RETURN 

1 FORMAT ! 9!1X, F5.2, F2..0), IX, F2.1, F3. 0/9! IX, F5.2,F2.0)/ 

1 3!1X,F5.2,F2.0) ) 

END 



RO 00010 
RO 00020 
RO 00030 
RO 00040 
RO 00050 
RD 00060 
RO 00070 
RD 00080 
RD 00090 
RO 00100 
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SUBROUT INE 



L L S C F T 



SUBROUTINE LL SOFT ( NVF , RMS , E RM , MODE ) 

C0MM0N/TABLE/N1,N2,N3,N4,N5,N6,F( 25, 12,ll,l,l,l),X(25,b),U<2t ,6) 
DIMENSION I (6J,N<6) ,G(26), A(3J 

EQUIVALENCE (1(1), ID, (1(2), 12), (1(3), 15), (1(4), 14), (1(5), 

1 15) , ( I (6) , 16), (N(l ),N1) 

SUBROUTINE LLSQFT PICKS OFF THE APPROPRIATE VALJtS FROM THE ARRAY 



LLSOOOlO 

LLS00020 

LLS00030 

LLS00040 

LLS00050 

LLS00060 

LLS00070 





F TO BE FIT BY WLSQ. 


LLS00060 

LLSCCC90 




ERM = 0. 


LLSOOlOO 




RMS = 0. 


LLSOOllO 




11 = 0 


LLS00120 


101 


11 = 11 + 1 


LLS00130 




12 = 0 


LLS00140 


102 


12 = 12 + 1 


LLS00150 




13 = 0 


LLS00160 


103 


13 = 13 + 1 


LLS00170 




14 = 0 


LLS00130 


104 


14 = 14 + 1 


LLS00190 




15 = 0 


LLS00200 


105 


15 = 15 + 1 


LLS00210 




16 = 0 


LLS00220 


106 


16 = 16 + 1 


LLS00230 




M = N(NVF) 


LLS00240 




J = 0 


LLSC0250 


110 


J = J + 1 


LLS00260 




KNVF) = J 


LLS00270 




G( J) = F(I1,I2, 13, 14,15, 16) 


LLS00280 




IF(J.LT.M)GO TO 110 


LLS0C290 




CALL WLSC(M,X(1,NVF) ,W(1,NVF) ,G,A,R,E) 


LLS00300 




RMS = RMS + R 


LLS00310 




ERM = AMAX1(£RM,E) 


LLS00320 




IF(MCDE.EQ.O)GO TO 130 


LLS00330 




DO 115 J=1,M 


LLS00340 




KNVF) = J 


LLSC0350 


115 


F( 11,12,13, 14,15, 16) = G(J) 


LLSC0360 




DO 120 J=l,3 


LLS00370 




KNVF) = J + N(NVF) 


LLS00380 


120 


F(I1,I2,I3,I4,I5,I6) = A(J) 


LLS00390 




KNVF) = M 


LLS00400 


130 


IF(I6.LT.N6)GG TO 106 


LLS00410 




IF( I5.LT.N5)G0 TO 105 


LLSC0420 




IF( I4.LT.N4)G0 TO 104 


LLS00430 




IF(I3.LT.N3)G0 TO 103 


LLSC0440 




IF( I2.LT.N2)G0 TO 102 


LLSC0450 




IFdl.LT.NDGO TO 101 


LLS00460 




wTS = 0. 


LLS00470 




DO 140 J=1,M 


LLS00480 


140 


WTS = WTS + W(J ,NVF) 


LLS00490 




NTS = WTS*N1*N2*N3*N4*N5*N6/M 


LLS00500 




RMS = SCRT(RMSZNTS) 


LLS00510 




RETURN 


LLS00520 




END 


LLS00530 



FUNCTION PHI 

FUNCTION PHI(I,X) 

GO TO ( 100,200, 300) , I 
100 PHI = 1. 

RETURN 
200 PHI = X 
RETURN 

300 PHI = X*X 
RETURN 
END 



PHIOOOlO 

PHI00020 

PHI00030 

PHIC0040 

PHI00050 

PHI00060 

PHI00070 

PHI00G80 

PHI00090 
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SUBROUTINE WLSQ 

SUBROUTINE WLSQ ( M , X , rt* F , A , R , E J WLSOOOlO 

DIMENSION X(l) ,F( 1) ,W(i) ,A( 1) WLS00020 

WLS00030 

SUBROUTINE wLSO OBTAINS A WEIGHTED LEAST SQUARES APPROXIMATION OF WLSOOQAO 
THE FORM A(1)*PHI (1,X) + A ( 2 ) *PHI ( 2 , X ) + A ( 3 ) *PH I ( 3 , X ) WLS00050 

WLS00060 

THE OALOULATION IS PERFORMED BY FIRST OBTAINING AN ORTHOGONAL WLS00070 

BASIS PF1,PF2,PH3 AND THEN FINDING THE APPROXIMATION WITH RESPECTWLS00080 



TO THESE FUNCTIONS. FINALLY, THE 
FOUND. 

PHKXJ = PHI(1,X) 

PH2(X) = PHI(2,X) - Q1*PHI(1,X) 
PH3(X) = PHI(3,X) - Q2*PHI(2,X) - 
A(2) = 0. 

A(3) = 0. 

01 = 0 . 

D1 = 0. 

B = 0. 

DO 100 I=1,M 
FT = PHI (1,X( I) ) 

01 = 01 + W(I)*FT*=i‘2 
D1 = D1 + W ( I )*FT*F( I ) 

100 B = B + W( n#FT*PHI (2,X( I ) ) 

Ad) = 01/01 
Q1 = B/Cl 

02 = 0. 

D2 = 0. 

B2 = 0. 

B3 = 0. 

DC 120 1 = 1, M 
FT = PH2IX(I)) 

02 = 02 + W(I)*FT*#2 
D2 = D2 + W(I)*FT*F( I) 

FS = PHI (3,X(I) ) 

B2 = B2 + W(I)*FT*FS 
120 B3 = B3 + W( I )=f'FS^PHI(l,X( I) ) 
IF(C2. EQ.O. )G0 TO 300 
P2 = 02/02 
Q2 = B2/C2 
04 = B3/C1 - Q1*Q2 

03 = 0. 

D3 = 0. 

DO 140 1 = 1, M 
FT = PH3(XU)) 

03 = 03 + W(I)*FT**2 
140 C3 = 03 + W(I)*FT*F(I) 

IF(C3.EC.O. »G0 TO 210 

P3 = 03/03 

A(l) = A(l) - P3*Q4 

A(2) = -P3*C2 

A(3) = P3 

210 A(l) = A(l) - P2*Q1 

A(2) = A(2) + P2 

300 R = 0. 

E = 0. 

DO 310 1=1, M 
FV = 0. 

DO 305 J=l,3 

FV = FV + A(J)*PHI( J,X(I ) ) 

305 CONTINUE 

F( I) = F(I ) - FV 
£T = ABS(F( I ) ) 

R = R + ET**2 
310 E = AMAX1(E,ET) 

RETURN 

END 



COEFFICIENTS OF THE PHI'S ARE WLS00090 

WLSOOlOO 

WLSOOllO 

WLS00120 

WLS00130 

Q4*PHI(1,X) WLS00140 

WLS00150 

WLSC0160 

WLS00170 

WLS00180 

WLS00190 

WLS00200 

WLS00210 

WLS00220 

WLS00230 

WLS00240 

WLS00250 

WLSC0260 

WLS00270 

WLS00280 

WLS00290 

WLS00300 

WLS00310 

WLS00320 

WLS00330 

WLS00340 

WLS00350 

WLS00360 

WLS00370 

WLS00380 

WLS00390 

WLS00400 

WLS00410 

WLS00420 

WLS00430 

WLSC0440 

WLSC0450 

WLS00460 

WLS00470 

WLSC0480 

WLS00490 

WLS00500 

WLS00510 

WLS00520 

WLSC0530 

WLS00540 

WLS00550 

WLS00560 

WLS00570 

WLS00580 

WLSC0590 

WLS00600 

WLS00610 

WLS00620 

WLS00630 

WLSG0640 

WLS00650 

WLS00660 

WLS00670 
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A. 6 Sample Data and Output for LSTSQ 



Essentially the same data was run on LSTSQ as on MINMAX. Va] es for 
the abscissae are necessary, and were simply given as equally spaced. 

The version of RD given in the listing was modified to read fewer than 
10 entries in the first variable by replacing the read statement by 
100 READ(5,1) (F(I1, 12, 13, 1,1,1),X(I1,1) ,11=1,4) ,X(I2,2) ,X(I3,3) . 

The input data and output follow: 



(Input data for LSTSQ, as punched) 



4 




3 


2 


1 


1 


1 




1 


0 


5 






1 . 




1 . 






.5 




1.0 


C 


2.0 


1 


3.5 


2 


4.0 


3 


0.0 


0 


-.5 


0 


1.5 


1 


2.9 


2 


3.8 


3 


1.0 


0 


0.0 


0 


1.5 


1 


2.2 


2 


2.5 


3 


2.0 


0 


- 1 . 


0 


0.0 


1 


1.5 


2 


2.5 


3 


0.0 


1 


- 2 . 


0 


-.5 


1 


.05 


2 


2.0 


3 


1.0 


1 


- 1.5 


0 


-.2 


1 


0.3 


2 


1.9 


3 


2.0 


1 
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(Output for LSTSQ ) 



THE INDEPENDENT VARIABLES 



WTl 


5.000E- 


•01 


l.OOOE 


00 


l.OOOE 


00 


5.000E- 


-01 


XI 


0.0 




l.OOOE 


00 


2.000E 


00 


3.000E 


00 


WT2 


l.OOOE 


00 


l.OOOE 


00 


l.OOOE 


00 






X2 


C.O 




l.OOOE 


00 


2.000E 


00 






WT3 


l.OOOE 


00 


l.OOOE 


00 










X3 


C.O 




l.OOOE 


00 











THE DEPENDENT VARIABLEt FIRST INDEPENDENT VARIATION ACROSS THE PAGE 



2 


13 


14 


15 


16 


















1 


1 


1 


1 


1 


l.OOOE 


00 


2. OOOE 


00 


3.500E 


00 


4. OOOE 


00 


2 


1 


1 


1 


1 


-5.000E- 


-01 


1.500E 


00 


2.900E 


00 


3.800E 


00 


3 


1 


1 


1 


1 


0.0 




1.500E 


00 


2.200E 


00 


2.500E 


00 


1 


2 


1 


1 


1 


-1 .OOOE 


00 


0.0 




1.500E 


00 


2.500E 


00 


2 


2 


1 


1 


1 


-2.000E 


00 


-5.000E- 


01 


5.000E- 


•02 


2. OOOE 


00 


3 


2 


1 


1 


1 


-1.500E 


00 


-2. OOOE- 


01 


3. OOOE- 


•01 


1.900E 


00 



THE ERROR ARRAY FOLLOWS 



12 13 14 15 16 
11111 



2 

3 

1 

2 

3 



1 

1 

2 

2 

2 



1 

1 

1 

1 

1 



1 

1 

1 

1 

1 



1.364E-01 -2.045E-01 
-9.094E-03 
-3.637E-02 
9.091E-02 
-2.136E-01 
-1.727E-01 



2.045E-01 -1.364E-01 
1.364E-02 -1.363E-02 9.090E-03 

5.455E-02 -5.454E-02 3.636E-02 

9.091E-02 -1.364E-01 1.364E-01 -9.091E-02 

3.205E-01 -3.205E-01 2.136E-01 

2.591E-01 -2.591E-01 1.727E-01 



THE FIT ON THE ITH VARIABLE, 



WITH RMS 
MAXIMUM 



ERROR = 
ERROR = 



1.935E-01 

3.205E-01 



THE COEFFICIENTS 
ALL BUT THE ITH 
THROUGH THE 7TH 



12 13 14 15 
1111 



2 

3 

1 

2 

3 



1 

1 

2 

2 

2 



1 

1 

1 

1 

1 



OF PHIi, PHI2, AND PHI3 ARE FUNCTIONS OF 
VARIABLE, AND ARE LISTED HERE AS THE 5TH 
VALUES FOR II, RESPECTIVELY. 



16 

1 



1 

1 

1 

1 

1 



8.636E-01 
1 -4.909E-C1 
1 3.637E-02 
1 -1.091E 00 
1 -1.786E 00 
1 -1.327E CO 



1.466E 00 -1.250E-01 
2.252E 00 -2.750E-01 
1.709E 00 -3.000E-01 
1.227E 00 3.576E-07 
8.534E-01 1.125E-01 
7.932E-01 7.500E-02 
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The output is largely self explanatory. The input function and 
error array are output in a format similar to that of Tables 2-6. 

The coefficients of the desired approximation are stored in the same array 
as the input function. On output, when the approximation was in the 
j variable, the coefficients of (l)^* and (f)^ are output as the 

(N^ + 1)^^ , (Nj + 2)^^ and (N^ + 3)^^ values of the variables. 

This is explained in the output, also. 
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